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Abstract 

We consider triholomorphic maps from an almost hyper-Hermitian manifold 
into a (simply connected) hyperKahler manifold This notion entails 

that the map u G satisfies a quaternionic del-bar equation. We work un¬ 
der the assumption that u is locally strongly approximable in by smooth 

maps: then such maps are almost stationary harmonic, in a suitable sense (in the 
important special case that A4 is hyperKahler as well, then they are stationary 
harmonic). We show, by means of the bmo— /t^-duality, that in this more general 
situation the classical £-regularity result still holds and we establish the validity, 
for triholomorphic maps, of the W^d_conjecture (i.e. an a priori W^’^-estimate in 
terms of the energy). We then address compactness issues for a weakly converging 
sequence U£ u^o of strongly approximable triholomorphic maps : M ^ Af 
with uniformly bounded Dirichlet energies. The blow up analysis leads, as in the 
usual stationary setting, to the existence of a rectifiable blow-up set E of codi¬ 
mension 2, away from which the sequence converges strongly. The defect measure 
0(a:)’H^”^“^ LE encodes the loss of energy in the limit and we prove that for a.e. 
point on E the value of 0 is given by the sum of the energies of a (finite) num¬ 
ber of smooth non-constant holomorphic bubbles (here the holomorphicity is to 
be understood with respect to a complex structure on Af that depends on the 
chosen point on E). In the case that A4 is hyperKahler this quantization result 
was established by C. Y. Wang uni with a different proof; our arguments rely on 
Lorentz spaces estimates. By means of a calibration argument and a homological 
argument we further prove that whenever the restriction of En(AJ \Sing^^) to an 
open set is covered by a Lipschitz connected graph, then actually this portion of 
E is a smooth submanifold without boundary and it is pseudo-holomorphic for a 
(unique) almost complex structure on AA (with 0 constant on this portion); more¬ 
over the bubbles originating at points of such a smooth piece are all holomorphic 
for a common complex structure on Af. 


1 Introduction 

Triholomorphic maps have appeared in high-dimensional gauge theory and string 
theory [T^ and play an important role for the geometric understanding of hyperKahler 
manifolds m From a more analytic perspective, they constitute an important class 
of harmonic maps and their more peculiar structure may lead to deeper results than 
those coming from the general high-dimensional harmonic map theory. 

Consider two compact oriented Riemannian manifolds {AA , g) and [Af, h) respec¬ 
tively of dimensions 4m and 4n and assume that the manifold A4 [respectively A/] 
carries three smooth almost complex structures i, j, k compatible with the metric g 
[respectively I,J,K compatible with h] satisfying the quaternionic relation ijk = —Id 
[respectively IJK = —Id]. To each of these almost complex structures, take for ex¬ 
ample i, we can associate uniquely a non-degenerate two-form Wj on AA defined by 
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III tli6 same way we have ujj and on Ai and (using the metric h 
instead) O/, rij, Qk on AA. In the case when Ai and A/" are hyperKahler manifolds all 
the almost complex structures are actually parallel complex structures and all of the 
two-forms defined above are closed. We will require, for the purposes of this work, that 
the forms 11/, fij, Qk on AAare closed, whilst we will not need the closedness for coi, u)j 
and Wfc on A4. A manifold Ai of the type just described is usually called almost hyper- 
Hermitian. By a result of N. Hitchin m Lemma 6.8] the closedness assumption on 
the forms on M amounts to the apparently stronger condition that A/" is a hyperKahler 
manifold. We will always assume that Af is simply connectec0 

By the Nash embedding theorem we can isometrically embed Af in for some 
large enough Q, and from now on we will assume to have done this. We say that 
u : Ai —)• AA is W^'‘^{Ai,M) if u is of class as a map from Ai into and 

moreover u{x) G AAfor ?^^”*-a.e. x G Ai. We consider u (and later on a sequence {u/}) 
of class W^'‘^{A4,Af) satisfying the triholomorphi(0 map equation 

du = Idui + Jduj + Kduk, (1) 

i.e. for a.e. x € A4 and any tangent vector Ai at x it holds 

du{X) = Idu{iX) + Jdu{jX) + Kdu{kX), 

where I, J, K act on A/". The situation that we will address is the following. We 
assume that the map u satisfies the following condition 

d{u*a) = 0 for every closed 2-form a on J\f. (2) 

Condition Q is rather easily seen to be true in the important case that u is locally 
strongly approximable in by smooth maps. In view of [6j we have that ([^ is 

actually equivalent to the local strong approximability when is torsion-free. 

This is satisfied for example when n = 1, i.e. A/" is a K3 surface. We are not aware of 
any example with non-trivial torsion in H 2 {Af, Z). 

A map u as above turns out to be stationary harmonic (i.e. a critical point of 
the Dirichlet energy for local perturbations both in the target and in the domain, this 
notion was introduced by R. Schoen m) when the forms cjj, coj, ujk are closed, i.e. 
when Ai as well is hyperKahler. The assumption ([^ is crucial her^as we will see. In 
the more general case that we address (i.e. LOi, ooj, ojk are not necessarily closed) then 
we find a suitable notion of “almost stationarity”: for any vector field X G (7“ (Ad) 

^The assumption that the hyperKahler manifold M is simply connected is not very restrictive: 
indeed, when it is not so, the manifold admits a finite cover that is the product of several simply 
connected hyperKahler manifolds with a flat torus T"^". The flat torus does not admit any harmonic 
S'^’s so for the purposes of the present work, where we will deal with compactness issues, it represents 
a very easy case as it does not allow bubbles to form. 

^With our convention for the PDE it is clear that (i, —/)-pseudoholomorphic maps are special 
cases of triholomorphic maps. The same goes more generally for {(ai -\- bj + ck), —(al -\- b,J -\- cK))- 
pseudo holomorphic maps, with = 1. However, as noted with explicit examples in [7], the 

triholomorphic notion is stricly more general than the pseudoholomorphic one. 

^Actually we will need ([^ only for a = fl/, Qj, Q.k, and these three forms do not generate the whole 
of H’^(A,Z), however we are not sure whether this actually yields a weaker assumption or whether 
the fulfilment of ([2| by a triholomorphic map is implied by its validity on these three forms only. 

^This was overlooked in [7]. 
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and denoting with = Id + tX the local variation induced by X then we have 


dt 


t=o 


I\V{uo^,)\^>-C I \X\\Vu\\ 
M M 


where C = ||(ia;j||oo + ||dwj||oo + ||dwA:||oo- A similar statement holds for local pertur¬ 
bations in the target. We observe further that u satishes a perturbation of the usual 
harmonic map PDE, namely (denoting by A the second fundamental form of N) 


AgU + A{Xu, Xu) = /(x, u,Xu ), 


(3) 


where / is smooth and has linear dependence on the hrst derivatives of u (observe that 
this type of perturbation is sub-critical). Moreover we prove in Proposition!^ an almost 
monotonicity formula for the energy ratio ^ 4 ^,^ |VmP at all x G At. 

Inspired by the similarities with stationary harmonic maps, we show that the basic 
results that are known in that classic setting are still valid in our situation, namely 
we prove an e-regularity result (see Propositions]^ and Q, i.e. there exists a threshold 
eg > 0 depending only on the geometric data (and not on the particular u satisfying Q 
and @) such that whenever the energy ratio of u is below Sq on a ball Bji{x) then u 
is smooth on B^j 2 {x) with a sup-bound on the modulus of its gradient. Moreover the 
singular set of u, i.e. the closed set where u fails to be smooth, is a p^^^^^-negligeable 
set. These results are based on the PDE analysis carried out in Sections 3.1 and 
where we hnd that the PDE ^ is of the form AgU G (where denotes the 
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inhomogeneous Hardy space) with Hardy norm controlled by the Dirichlet energy of u. 
Then we use the bmo — duality, similarly to the case treated by L. C. Evans |in) . to 
infer the e-regularit}0 We also establish the -conjecture for triholomorphic maps, 
i.e. the analogue of the a priori VP^’^-estimate conjectured to hold by T. Riviere 
page 5] in the case of stationary harmonic maps: 

Theorem 1.1. Let u : Ai —)• W C 


be a triholomorphic map, i.e. (Til) holds. Then 


u G W'^'^{M) and \\u\\y/ 2 ,i(^j^^ < Consti^j^^g-^ {Const(^j^^^Q) -\-E{u)^ , 


where E{u) denotes the Dirichlet energy of u and Const^^jy^-^g-^ is a constant depending 
on the embedding of M and Constf^j^ g-^ depends only on the geometry of Ai. 


rem 


This will be a fundamental ingredient in the proof of our quantization result Theo- 
1.2 In the case of stationary harmonic maps from a domain in 


into a spherical 

target the IT^’^-estimate is known to hold thanks to the observation by F. Helein that 
Au G the homogeneous Hardy space (see |9], |19jl. This a priori estimate is the 
starting point of the proof in | 10 | . where the BMO — iA^-duality is exploited. 


After these preliminary results we proceed further with compactness questions for 
a sequence {u^} of triholomorphic maps (with assumptions Q and Q as above for 
each Ui) with uniformly bounded Dirichlet energies. Without loss of generality we 
may assume Ui ^ u weakly in W^’^. It is natural to ask what knowledge we have 
of the limit map u and of its smoothness properties and if (and on which set) we 

®Evans proved the result for stationary harmonic maps when the target is a sphere. For stationary 
harmonic maps with arbitrary target the £-regularity result was established by F. Bethuel [5]. For 
minimizing harmonic maps it had been earlier proved by R. Schoen and K. Uhlenbeck I33- 
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may expect a stronger form of convergence. The qnestion is wide and open to a large 
extent (even in the case where both M and M are hyperKahler) and its solntion wonld 
immediately impact the compactification of the moduli space of triholomorphic maps, 
leading fnrther to the dehnition of geometric invariants of hyperKahler manifolds by 
means of a snitable connting of triholomorphic maps, as ontlined by the second anthor 
in m- The identihcation of a snitable class of triholomorphic maps is one of the tasks 
here, and it is also in this light that condition Q shonld be viewed. In Section we 
show that the fnndamental resnlts of F.-H. Lin |23j for stationary harmonic maps hold 
with corresponding statements in onr case, namely: we can identify a closed set S 
of locally hnite ^^^"^“^-measnre away from which the convergence {n^} —)■ n is strong 
(even C^), the set S is (4m — 2)-rectihable, the Radon measnres |weakly- 
* converge to the measnre \_T, with density Q > Sq a.e|^ 

Moreover u is smooth away from a closed set Sing^ C S (on which it is discontinnons) 
and satishes globally the PDF A.gU + A{Vu,Vu) = f{x,u,'Vu). However, whilst u 
is smooth away from Sing„ (and thns u is almost stationary on A4 \Sing„) the same 
might a priori not hold on the whole of Af. A priori we conld have that (for some 
closed 2-form a) the form d{u*a) is a non-zero distribntion with snpport on the set 
Sing^, since the weak VF^’^-convergence does not allow to pass condition ([^ to the 
limit. The nnderstanding of whether the local strong approximability holds across 
Sing^ is of central importanc^ It is very mnch related to the nnderstanding of the 
strnctnre of S as well. To see this, let ns focns on the important case where A4 and 
N are hyperKahler: the lack of the strong approximability condition wonld imply a 
lack of stationarity of u and (by the blow np formnla in [22]) a lack of stationarity for 
the varifold S, localized on Sing„, with obvions implications on the regnlarity. Remark 
that a priori S might as well have bonndary contained in Sing^. 

The main compactness resnlts we provide are the following two theorems. The hrst 
regards the quantization of energy for the defect measnre: 

Theorem 1.2. Consider a sequence of maps ui : A4 ^ M (with M almost hyper- 
Hermitian and M hyperKahler of arbitrary dimensions respectively 4m and 4n ) satisfy¬ 
ing and with uniformly bounded Dirichlet energies. Let ui ^ u weakly in W^’‘^ 
and let S be the blow-up set, i.e. ^ \'S/u\‘^d'H'^^ + Q{x)d'H^'^~'^ CT, with 

0 > Eq > 0. Then for 'H^^~‘^-a.e. x ^ T, we have 

e{x) = f2E{(f,), (4) 

where G N and each (fg : S'^ ^ Af is a smooth non-constant harmonic map (these 
are called bubbles). 

In other words we prove that the whole loss of energy encoded in the defect measnre 
comes from bnbbling off of S'^’s and there is no leftover energy in the necks connecting 
the bnbble^ For L^'^™'“^-a.e. x G S we fnrther have that the babbles are holomorphic 

®The measure L E is usually called defect measure and encodes the possible loss of 

energy in the limit and lack of strong convergence. The set E is usually called blow-up set, bubbling set 
or concentration set. It is the set where the energy concentrates, leading to lack of strong convergence 
and formation of bubbles. 

^Unfortunately [7] does not notice that it plays a role In the analysis of the regularity of E. 

®A quantization result of this type is not known for arbitrary stationary harmonic maps, in |24) 
the target is a sphere. 
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S'^’s for a certain choice of almost complex structures jx and Jx (depending on the 
point x!), see Proposition]^ This quantization result was proved, in the case where 
both M and A/" are hyperKahler, by C.-Y. Wang m- In the present work we extend 
that result to our more general geometric setting, using a different proof inspired by 
the techniques of Lin-Riviere ||24I |. In Remark 1.5 and Section]^ we describe a gauge 
theoretical application where the more general setting is required. 


Finer properties blow-up sets are generally a very difficult task and not much is 
known beyond rectihability. With the second compactness result we start, in the case 
of interest, a more detailed analysis of the structure of the blow-up set T, and provide 
a hrst important step towards the understanding of it. Roughly speaking we show, by 
means of a new idea based on a homological argument and a calibration argument (and 
using the quantization result), that whenever S \ Sing.^ has a bit of regularity, then it 
is pseudoholomorphic for a hxed almost complex structure and the bubbles produced 
are holomorphic spheres for a hxed complex structure in Af. Precisely we have: 


Theorem 1.3. Under the same hypotheses of Theorem 1.2. whenever we have an open 
ball C Ai \Sing^ such that S n is contained in a (connected) boundaryless 

Lipschitz (4m — 2)-submanifol^C C and ^'^”^“^(11 n > 0 then there exist 

constants a, b, c with -\-b'^ -\-(? = 1 such that 


(i) coincides with the whole of C and it is pseudo holomorphic with respect to 
the almost complex structure ai + bj -\-ck (a posteriori this pseudo holomorphicity 
implies that C was actually smoooth); 

(ii) for every x G S n B^'^ the bubbles ((/>*)* (5^) that appear in ^ are holomorphic 

with respect to the complex structure —{al -\-hJ -\- cK) on Af. Moreover 0(x) in 
is constant on S n 


Remark 1.1. The result of |22| gives that, in the case when A4 is also hyperKahler (and 
thus the maps ue are stationary harmonic), the varifold 0(x)(ip^^™'“^LS is stationary 
away from Sing^: well-known GMT results (see Allard’s theorem lU) yield that for a 
stationary varifold there exists a dense open subset where the varifold is smooth (but 
stronger regularity properties for arbitrary stationary varifolds are unknown, although 
expected). In particular, whenever u has a continuity point on S, the assumption on 
S in Theorem 1.3 is somewhere satished. The latter theorem thus yields that each 
connected smooth piece of S is actually (pseudo)-holomorphic (of codimension 2 in 
AA). The point here is that, for each such connected component, we have a hxed almost 
complex structure: this also implies that for all points in such a connected component 
the almost complex structure for which the bubbles are holomorphic is independent 
of the point (compare the weaker statement mentioned after the quantization result). 
To our knowledge this is the hrst instance of application of a quantization result for 
stationary harmonic maps in a high-dimensional situation. 


We cannot prove that the whole of S needs to be (pseudo)-holomorphic for a unique 
almost complex structure. A priori it might be conceivable, for example, that several 
smooth pieces (each one pseudo holomorphic for a different structure) come together 
along a submanifold of codimension 3 that is contained in Sing^. Ruling out such a 


®Here we mean that for every point on the submanifold there exists an open cylinder in which the 
submanifold can be expressed as the graph of a Lipschitz function from —>■ (upon rotating 

and relabelling the coordinate axes). 
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behaviour would of course impact the knowledge of the regularity properties of u. This 
is an aspect of the deep connections between the stationarity properties of u and the 
regularity properties of u and S. We conjecture the following. 

Conjecture: With the same assumptions as in Theorems 11.21 and 1.3 Sing„ has 


codimension at least 3 and S is made of a hnite number of (pseudo)-holomorphic vari¬ 
eties (each connected component is holomorphic for a certain almost complex structure 
on A4, possibly all of these structures might be the same), each of which is a cycle in 
Ai \Sing^ and whose common boundary is the (4m — 3)-dimensional stratum of Sing.^. 


Remark 1.2. A very important problem to address is the existence of a homogeneous 
triholomorphic map with a singular set of codimension 3. The only known example 
is constructed in the case of a non-compact target (this is shown in |8j based on the 
construction of 0 ). The non-existence of such maps would impact the previous conjec¬ 
ture, which could be strenghtned considerably: all varieties would be cycles and Sing^ 
would have codimension at least 4. 


Remark 1.3. Another important question to address, on which however we do not 
focus at all in this work, is the optimal estimate on the size of the singular set for a 
triholomorphic map satisfying Q: is it possible to show that the codimension is at 
least 4? We wish to stress here that the example of triholomorphic map with singular 
set of codimension 3 mentioned in the previous remark is stationary harmonic but does 
not satisfy ([^. 

Remark 1.4. We believe that the argument given in Section |^for the proof of Theorem 

for the moment we prefer to limit ourself to the given statement and postpone any 
improvement to future work. 

Remark 1.5. As we described above, we extend the triholomorphic notion to the case 
where A4 is almost hyper-Hermitian (rather than hyperKahler). We then show that the 
hrst order PDE Q and the geometric structures on the manifolds A4 and AAlead to a 
jacobian structure for AgU which enables the use of Hardy space techniques (this plays 
an important role both in the e-regularity result and in the quantization result). Relax¬ 
ing the assumptions on Ai is interesting not only as an example of “almost stationary 
harmonic maps”, as we described earlier. Indeed the analytic and geometric techniques 
employed in the proofs of our results can hnd a rather direct application in the problem 
treated e.g. by T. Walpuski |39| . In that work the author deals with triholomorphic 
sections in a bundle of hyperKahler manifolds (with m = 1, using our notations); the 
problem originates in Gauge Theory on Spin{7) manifolds. The PDE satished by these 
triholomorphic sections corresponds to our equation Q with in addition lower order 
perturbation terms (depending on u but not on its derivatives). Our analysis can be 
carried out in the same fashion in the application under consideration, in particular 
it gives an affirmative answer to the quantization question, which is conjectured to be 
true but left open in |391 page 5]. We postpone to the hnal section a brief explanation 
of the modihcations involved. 


1.3 is robust enough to be pushed to a situation with weaker assumptions, however 


Acknowledgements: The hrst author wishes to thank Tristan Riviere for very 
fruitful conversations while this work was in progress. The authors wish to thank 
Thomas Walpuski for bringing their attention to the topic of Fueter sections, described 
in Section [3 
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2 Dirichlet energy and almost-monotonicity formula 

We begin this section by recalling (see [7j, |23) that the Dirichlet energy density 
for a triholomorphic map can be written, thanks to ([^, as 

\\Vu\l = A u*ni + 1 A u*^j + ^ A u*nK] • 

(5) 

This holds pointwise a.e. on the domain for u G and it snfhces to prove 

it at a point x by choosing coordinates so that {TxM,i,j,k) is isometrically iden- 
tihed with where {_/df=i the nsnal qnaternion nnits and 

are the standard associated Kahler forms. Consider a standard orthinormal basis 
{ei,jiei,j 2 ei,j 3 ei, ■■■, The differential du is then a linear map into 

Tu[^)M satisfying du = Iduji + Jduj 2 + Kduj 3 . The compntation in [7j Proposition 
2.2 shows that for any map v, not necessarily triholomorphic, it holds 

- .o ^ A v*nj + (aa)'™-' A v*nj + (as)"’""' A v^Qk] = 

[2m — Ij! ■' 

= 2 ~ ~ ~ ~ ( 6 ) 
and in particnlar, for u triholomorphic, we hnd 

-Tr^-iAT A u*nj + (aa)'™-' A u*Qj + A = 

[2m — Ij! ■' 

= \\M^^ 

The signihcance of ([^ and ([^ relies in the following fact. Under the assnmptions 

(a) At and A/" are Hyperkahler (so a;,, cjj, a;^, D/, Dj, are closed) 

(b) whenever a is a closed 2-form on AAthen d{u*a) = 0, cf. (|^ 

the qnantity - (2m-i)! / A A u*klj + (0;^)^™“^ A is 

M 

a null Lagrangian for local variations in the target and in the domain (see e.g. |30) 
for the analogons statement and proof in the case of psendoholomorphic maps). Then 
nnder snch variations, for u triholomorphic, in view of ([^ we can see that u is weakly 
harmonic and stationary harmonic. In the case that u is a smooth triholomorphic map 
it follows from (|^ that u minimizes the Dirichlet energy in its homotopy class. 

We will be interested in a more general setting, namely we will drop the assnmptions 
that uii, ujj, ujk are closed. Then the same argnment shows that a triholomorphic map 
u is almost stationary harmonic, in the sense specihed in ([^ below. 

Let X G (7“ be a vector held on Ai and 'I't = Id + tX be the corresponding 
1-parameter family of domain variations. We know that 

y’|V(novI/i)|2> (7) 

M 
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^ A {u O ^tynj + A (n o ^tT^j + ^ A {u o ^1/*)*^!^] , 

[Zm — 1 ): '- ■' 

by (1|), with equality for t = 0 since u is triholomorphic. Now 

— (no = — 'I'it6*n7 = Cxu*Qi = d{ixu*^i), 

where we used Q in the last equality. Using Stokes theorem 

^ I A (no = I diiuif^-y A {ixu*ni). 

M M 

Doing the same for the remaining two terms we get 

^ J [(cn*)2—1 A {u o ^tyni + 1 A {u o ^tyQj + A {u o ^ty^x] = 

M 


= j d{{uj,y^-y A iixu*nj) + diicojf^-y a iixu*nj) + dHu^kf^-y a iixu*nK) 

M 

which is bounded in modulus by 


(2m - 1)! / C\X\\Vu\'^, where C = ||da;i||oo + ||dwj||oo + \\diOk\\oo- 


M 


Putting together the latter equation and ([^ we hnd 


d 

dt 


t=o 


j |V(noxI/,)|2 >-C J \X\\Vu\y 

M M 


( 8 ) 


If cjj, ojj, Uk are closed then (7 = 0 and we have that u is stationary harmonic. If (7 7 ^ 0 
then we have a bound on how much the energy can decrease. In other words Q gives a 
notion of almost-stationarity. Remark that by zooming around any point of Ai at very 
small scales, we can make (7 as small as we wish. This notion of almost-stationarity 
has the same flavour as the notions of almost area minimizing current or semicalibrated 
current in the setting of Plateau’s problem. 

Remark 2.1. We do not give the argument to show the weak harmonicity of u for 
local variations in the target Af when A4 is hyperKahler. The argument is 

similar to the above and does not even require assumption Q, compare also |30| . 
Moreover by dropping the closedness assumption on Ui, Uj, ujk we hnd once again almost 
stationarity w.r.t. local variations in the target, namely: for every (j) G 
and denoting with the nearest point projection from onto Af (well dehned in a 
tubular neighbourhood of A/) we have f |V (IIt/o (u + t(j))) p > —C f |</>||Vnp. 

M M 

Monotonicity properties. We will see next that the Dirichlet energy of a tri¬ 
holomorphic map u : Ai ^ Af satisfying ([^ has good almost-monotonicity properties. 
In the special important case that Ai is hyperKahler then u is actually a stationary 
harmonic: in this case the almost-monotonicity formula for the energy ratio of triholo¬ 
morphic maps follows from the corresponding result for stationary harmonic maps [28], 
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|33) . We present next a direct proof of this monotonicity formula, based directly on 
the first order PDE Q and that has the advantage of extending to the case that Jv[ is 
merely almost hyper-Hermitian. 

Proposition 1. There exists ro > 0 such that for any x € A4 and r < xq the following 
quantity is monotonically decreasing (weakly) as r f 0 ; 


(1 + (4m — 2)r) 


-4m—2 


! Br{x) 


A u*Tli + (as A u*Tlj + (a^™”^) A u*nK 


(9) 


This immediately translates into an almost-monotonicity formula for the energy ratio: 
there exists tq > 0 such that for r < xq and x € M 

= /(^) + O (r/(r)) with f a non-decreasing function of r. 

^ JBr{x) 

Proof. Case 1: Let us consider first the case when A4 is replaced by the unit ball in 
endowed with the flat Euclidean metric and with the standard complex structures 
71 ) 72)73 and associated Kahler forms oi, as, 03 . In other words the domain is the unit 
ball in H™, the m-dimensional quaternionic vector space. 

We will denote by xdx the radial vector field X]a=i ^adxai by dr the radial field of 
unit length for x = (xi,..., x^rn) ^ ( 0 , •••, 0 ), and by dr = d\x\ the differential of the 
function x (that is the metric dual to dr for the flat metric). Computing explicitly the 
Lie derivative, for each i G {1,2,3}, we get Cxd^cii = 2a£ and therefore, distributing 
the derivative on the wedge product, 

^xa.(a|—i) = (4m-2)(a^-i). 

By the Cartan formula £. = id dt we can see that the 1-form (a^™'”^) satisfies 


Denote by 
1 

(2m — 1)! 


d{:'xdx 


2m-i)) ^ ( 4 m_ 2 )(a 2 ™-i) for each £g{ 1,2,3}. 
the tangential part of that is (for \x\ 7 ^ 0 ) 



1 )! 


a?”-* 


dr A LQr 


1 

(2m- 1)! 


(a 


2m —1 

e 


— '^{dr A jidr')^ 


( 10 ) 

where * denotes the Hodge star operation combined with the metric duality vec¬ 
tors/covectors. For the sequel let us remark that (here |x| 7 ^ 0) 


d 


lo-|4m—2 


d(^xa.(a^-^)) 

|^|4m—2 


(4m - 2)d|x| A i,^a^(a^’^ ^) 

I j.|4m—1 


= (4m 


- d\x\ 


2 ) 


2m —1 
£ 


|4m-2 


(4m 


2) 


(a 


2m— 

i _ / tan 

j.|4m—2 


We can now start the computations for the monotonicity formula, denoting by H/j 
and Bs two balls cenetered at the origin with radii respectively R and s (with R > s) 
and keeping in mind that u*kli and ag are closed forms: 


1 

j^Am—2 



2m—1 
(«i 


A u*Tlj 


1 

gAm—2 



(a?™-! 


A u*TLi = 
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Iqb^ (4m - 2)i?4m-2 


Au*ni - 


A u*ni - 


I (o^l ) 

Iqb^ (4m - 2)s'^™-2 

r . /^2m—1\ 

f t'xdx (®1 ) 


IoBr (4m - 2)|x|4™-2 (4m - 2)|x|4™-2 

_ f 


= 

Au*^?/ = 


ld{BR\Bs) (4m-2)|x|4™-2 


A?x*J2/ = 




P / 2m—1\ 

I V 1 _ /tan A ^.*Q 

,B,\B/V(4m-2)|xr-2y''““--^^^^^ |xr-2 ^■ 

In the last step we used ([^. Observe now, in view of (10), that 

(2m - 1 )! < A i.T!,)(dvol"'") = (»-fi,)(S, AJA) = 


( 11 ) 


= 9.i{du{dr) A du(/’i9r-)) = {Idu{dr),du{jidr)h\^,,) = -{du{dr),Idu{jidr)h\^,^y (12) 


where h|„( 2 ,) denotes the metric in the target at the point u{x). 

Writing the analogues of (llip and (12) for the remaining two terms, i.e. those 


involving the structures (j2,J) and (j^,K), and adding up the three parts to get the 
Dirichlet energy, as indicated in ([^, we get 


1 


^4m-2 


|Vn|2- 


1 


^4m—2 


|Vn|2 = 


{du{dr), I du(jidr) + Jdu(j 2 dr) + Kduijsdr))) 


\x 


4m-2 


Ba 


Bs 


Br\Bs 


(13) 

Using the triholomorphic equation du = Iduji + Jduj2 + Kduj^ in the r.h.s. of ( |13[ ) we 
finally conclude 


JiL, I IV..P 

Br 


1 


o4m—2 


|Vnp = 


\du{dr)\‘^ 

\rp\im—2 


Bs 


Br\Bs 


(14) 


This is indeed the usual monotonicity formula for stationary harmonic maps having 
as a domain an Euclidean ball |33| . Of course there was nothing special in the choice 
made of centering the balls Bb and Bg at the origin, we can write the same formula 
for arbitrary centers, by replacing the radial vector field xdx emanating from the origin 
with one emanating from the new center. We remark that the classical method for 
obtaining the monotonicity formula relies on the stationarity of the map with respect 
to radial variations in the domain, see [33] ■ We have replaced this by the Lie derivative 
computation with respect to xdx, in order to exploit the first oder information (|^. The 
proof just given is in the style of the proof of the monotonicity formula for calibrated 
cycles [T8] . 


Case 2. In the following we will modify this argument to obtain an almost mono¬ 
tonicity formula for A4 almost hyper-Hermitian. Let xq € Ai and consider a geodesic 
ball around xq, taking normal coordinates so that we work (up to a dilation) in the 


10 


















unit ball of centered at the origin. By suitably choosing coordinates we can 
make sure that Wi(0), aJj(O), a;fc(0) coincide with ai, 02 , as. In this way ai is the 
parallel extension of a;i(0) to the unit ball (parallel with respect to the flat metric), 
and similarly 02 and as are the parallel extension respectively of Wj(0), Wfc(O). Since 
— ai|(x) < IIVa;j||oo|a:| we can assume, up to dilating enough, that on the unit ball 
Bi(0) we have 

\oJi — Q!i|(x) < e|x|, \ujj — a 2 \{x) < e|x|, \ujk — asK®) < e|ic| 
for some small e > 0 and similarly 

\i - ji\{x) < e\x\, \j - j 2 \{x) < e\x\, \k -js\{x) < e\x\ 

while for the metric, by the choice of normal coordinates, 

\g- go\{x) < e|xp, 


where qq is the Euclidean metric. In other words we are dealing with a small 
bation of the situation that we treated earlier. We choose e so that e < 
unversal constant C that will appear in (18). 


pertur- 
for the 


The hrst observation is that, since the Dirichlet energy is given by ([^, in view of 
the smallness assumptions just made, the quantity 


A u*ni + {a2f^-^ A + {asf^-^ A u*nK (15) 

is a good approximation of the energy density |Vn|g(x), i.e. their difference is bounded 
in modulus (pointwise a.e.) by Ce\x\\S/u\^{x) for some dimensional constant C. Simi¬ 
larly the quantity 


-1 




'B. 


(ai)"™-' A u*ni + (as)"'"'^ A + (as)"™"' A u*nK (16) 


is a good approximation of the energy ratio l^^lg > be. their difference is 

bounded in modulus by l^^lg higher order term in R). So the energy 

ratio and (16) have the same asymptotic behaviour, therefore we will aim for an almost 
monotonicity formula for (16). This will be achieved through a modihcation of the 
argument given in the flat case above. 

A key observation for this modihcation is to observe that, in the present situation, 
although the following expression, i.e. the r.h.s. of (|11[), 


A “-fi/ + A A u-Qk (17) 

is no longer > 0, we still have a good lower bound, thanks to the fact that we are 
handling small perturbations of the previous case. Let us see why. As before we have 
that this expression is 


(expression (17)) 


{du{dr), Iduijidr) + Jdu(j2dr) + Kduijsdr)) = 


= {du{dr), Idu{idr)+Jdu{jdr)+Kdu{kdr)) — {du{dr), Idu{i—ji)dr+Jdu{j—j2)dr+Kdu{k—js)dr) 


= {du{dr), du{dr)) — {du{dr), I du{i — jl)dr + Jdu{j — j2)dr + Kdu{k — js)dr) > 
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(18) 


> \du{dr)\^ - Ce|x||Vu|^, 

for some universal constant C. This is the desired lower bouncf^ 
In order to complete the proof we remark that 


\ _ - 3)d|x| A 


, 2m—1 


2m—1\ 


d 


|4m—3 


14m—3 




4m—2 


(4m —2)(a^'” [Am — 2>)d\x\/\ iQ^ag, ^ + (4m — 3) (a^”* 

.14771 — 3 


tan 


14m—3 


\X\ 


14m—3 


= \X\ 


^2m-l ^^_3)(a“) 


tan 


|4m—2 


|4m-2 


(19) 


This suggest^^that a suitable choice for the “r.h.s. of a monotonicity formula” is (here 

|x| / 0) 


/ 2m—1\ / 2m—1\' 

* (4m - 2)|x|4^-2 + ^ ^ = 


(l + (4m-3)|x|)(a?”^ ^)^^^ \x\al 


2m—1 


|4m—2 


+ 


|4m—2 


An*H/ + (a2, J)-term + (Q;3, -fC)-term > 


> 


(1 + (4m — 3)|x|) [\du{dr)\‘^ — Ce\x\\Vu\‘^g) + |x|(l — CelxDlVul 


X 


4m—2 


> 


> 


(1 - Ce) \x\\Vu\l - (^(4m - 3)C + Cj e|x|2|Vn| 

I ^|4m—2 


> 0 


for e < and |x| < rg for a small enough rg. 

We can now integrate by parts as follows (here Bg and Bji are balls centered at 0 
of radii s < R < tq) 


0 < I d I ... J _|_ —^ \ ) A u*Qi + (0:2, J)-term + (03,il')-term 


Br\Bs 


(4m —2)|x|^™' ^ ^ 




/ 2m—1\ / 2m—1\' 

(4m - 2)V|4>n-2 + + («3, i^)-term = 


(1 + (4m — 2)R) 
(4m - 2)i?4m-2 


(«! ) A u*Qi + (02, J)-term + (03, iir)-term 




^In view of this lower bound it is natural to attempt to replace expression (171 in the r.h.s. of (131 
with something of the type ’expression 0 + .^1*1 expression (|15[)’, which is non-negative for K large 
enough. However we want, In addition, that this non-negative quantity should allow an “integration 


by parts- trick” as in (111. For this reason we will need to modify this tentative quantity slightly. 
^^This is how we modify the tentative quantity from the previous footnote. 
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(1 + (4m — 2)s) 
(4m — 2 ) 5 ^™-“^ 



^ u*^l + {o( 2 , J)-term + ( 03 ,/sr)-term 


I A u*Qi + {a 2 , J)-term + ( 03 , i^)-term — 

Br 

I (a^™“^) A u*^j + {a 2 , J)-term + ( 03 , /C)-term . 

Br 

This provides the desired monotonicity formnla (|^ , from which the almost monotonic¬ 
ity of Proposition follows immediately. 

□ 

Immediate consequences of almost monotonicity. The almost monotonicity 
formnla implies that for every point x € M there is a well-defined limit of the energy 
ratio as r —)• 0. This is called energy density of n at x and denoted 

by Q(u,x). A standard argnment exploiting the almost monotonicity also gives that 
Q{u,x) is an npper semi-continnons fnnction of x. 

Moreover there exist rg > 0 and C > 0, depending only on the geometric data, snch 
that whenever x £ M. and r < ro we have > &{u,x) — Cr. 

3 Analysis of the PDE 

3.1 Jacobian structure 

Proposition 2. Let u : A4 —)• AA C he a triholomorphic map in 

and Af Then AgU = {f^,f^), where the f^’s are in the Hardy space 

and ||/|Ui(w) < C/lVup for U C A4, where C > 0 is a univeresal constant (i.e. 
U 

independent of the particular u). Here Ag denotes the Laplace-Beltrami operator on 

{M,g). 

proof of Proposition It suffices to prove the result on a chart U <Z Ai. For any 
y £ Af C we can extend the almost complex structure I (dehned as a linear 
automorphism on TyAf) to a linear automorphism Igxt of TyMP = For that it 

suffices to decompose = Ty Af ®{Ty Af)'^ and define the new automorphism Axt 

to agree with I on whilst Lf>xt = 0 on {TyAf)'^. In the same way we can extend 
J and K to Jext and Kf.xt respectively. Introducing local coordinates x^ on Ai for 
I £ {1, ...4m} and standard coordinates ?/“ on for a £ {1, ...Q}. With these 
coordinates in mind the almost complex structures on A4 and the extended ones in the 
target will be expressed through coefficients as follows (we write it for j and Text only, 
the same notation will be used for i, k, Axt, Kext) 

^ {r.-Q}. (20) 

will be using greek indexes on the target and latin indexes on the domain, avoiding in the 
latter the use of the indexes i,j, k since these letters denote the almost complex structures on M. 


(1 -h (4m - 2)R) 
(1 A (4m — 2)s) 

^4771 —2 
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with implicit summation over repeated indexes. Remark that the extended automor¬ 
phisms are not almost complex structures in TyMP as they are degenerate in the direc¬ 
tions normal to T M. 

Since on the r.h.s. of Q we have I, J, K acting on vectors that are tangent to Af 
we can substitute lext, Jext, K^xt- Therefore the map u, as a map u = ...,v9) into 

satishes the equation 


du = lextdui + Jextduj + Kextduk. 


We will now rewrite this equation in the chosen coordinates. Using (20) we can rewrite 

f d \ 

Jextduj i \ = J^ji ■ 


,du'^ d 


dx^ dy^ 

and from here we have the triholomorphic map equation in the form 
du^ 






( 21 ) 


dx^ dx^ ' 

We will denote, as it is customary, by Qab the metric tensor g in local coordinates 
and by its inverse. The Beltrami-Laplace operator div{W) on {AA,g) is expressed 
in local coordinates by (always summing over repeated indexes) 


A5 = 




where ^/\^\dx^ A ... A dx^"^ is the volume form on AA. Applying this operator to 
(for a hxed fd G {1, ...Q}) and using (21) we hnd 


1 


Vgdx' 




^ _a 

^\dx' 

1 d 




dxS Je 


( 22 ) 


+ 




dx^dx^ 


The summation over 1=1, ...,4m in the last expression only involves products of the 
type g°‘^ig, g°‘^k^. Recalling the dehnition of cjj = g~^i (the non-degenerate two- 

form uniquely associated to the metric g and the almost complex structure i) we can 
see that for each (a, s) the term g°'^i^ is the coefficient of coi corresponding to dx°‘ ® dx^ 
when we express oji in the local coordinates, i.e. = oji ® skew- 

symmetry of u!i implies that g°‘^i^ = The same goes for the sum of products 

AS ^aius 

y y 


m 


Using the skew-symmetry w.r.t. (a, s) of g°'^i\, g°‘^j£, g°'^k^ and the symmetry i 

ct 

{a,s) of we can see that, summing over all (a,s), the second, fourth and sixth 

terms of (22) vanish. Using the skew symmetry of g^'^i^, g^'^j'i, g°‘^k^ again we can 
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rewrite the first, third and fifth terms of of (22) summing on (a, s) with s < a only. 
These observations lead to 




P - 


\/H 


E 

s<a 


d 

dx^- 


s\a“ii 






+ 


V\9\ 

1 

\/W\ 


The coefficients /( 


s<a 

E 

s<a 

p 


_d 

dx' 


du°‘ 

dx^ 

clu“ 

dx^ 

du°^ 

dx^ 


_ 

dx^ 

d 

dx^ 

JL 

dx^ 




dx°‘ 




du'^ 


(vIsI/'tfVd) 


(9x“ 

du°‘ 


+ 


+ 


(23) 




a , 


Ja, Ka are evaluated at u{x). The functions /( 


0 
a, 


J0 

5 


Ka are 


, so they can be extended 
- as smooth functions. As 


smooth on M (which is a compact closed submanifold of 
to a tubular neighbourhood first and then to the whole of 
u G lT^’^(Ad,M*^), the functions Ia°u, Ja°u, Kaou, that appear in the r.h.s. of (23), 
are therefore in lT^’^(Ad, M). Setting (for each (a, s)) ijj = ^/\g\g°'^ifIa and (j) = we 
can recognise, in the r.h.s. of (23), a sum of terms of the type , for ^ and 

4) in lT^’^(Ad,M). The Jacobian structure of this expression is the crucial assumption 

in the Hardy space /t^ (Af). Moreover the 


to infer, from [9], that lies 

Hard^orm of is control 
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ec 


by ||ViA||i 2 ||V</-|U 2 . So we can conclude 
that (23) is a PDE of the form KgU^ = / £ /t^(AJ) with ||/||/;i < Const(^j^ f^E{u). 

□ 


A classical result, e.g. |36| |31| . gives that any solution to An = / with / in the 
Hardy space /t^ belongs to and ||n||^y 2 .i < Const\\f\\fii. In the next subsection 

we show that the same holds for our map n by adapting the result for the classical 
Laplacian to the Laplace-Beltrami operator by a perturbation argument. We will thus 
establish the IT^’^-estimate of Theorem ll.il 


3.2 Hardy space and proof of Theorem |1.1| 

Let $ be a function in the Schwartz class 5(M^™') of rapidly decaying smooth 
functions such that J = 1. For t > 0 denote by the function The 

inhomogeneous Hardy space /t^ defined in terms of a maximal function: given / G 
^i(]^4m) classical Hardy-Littlewood maximal function is Mf{x) = sup^ |/b /|> 
with the supremum taken over all balls centered at x. We will consider a modification 
of this, namely supg^^^^ \ (roughly speaking the balls in the classical maximal 

function cannot get as large as we want - this is a local version of the Hardy-Littlewood 
maximal function). The inhomogeneous Hardy space is the following subset of 
Li(M^™) 


/fi = {/ E : sup * /I G L^M^™)} (24) 

0<t<l 

^^One of the versions of the result [9], that applies to our situation, is the following. We state it in 
R™, taking it from |20) . Section 4.2. Let a € 1T^’^(R’") and /3 be a closed (in distributional sense) 
(m — l)-form of R"* with coefficients in L^fR"*); then da A j3 = fdx^ A ... A da:"*, where / belongs 
to the Hardy space "H^fR"*). From this we see the following: given af) and 4> in 1T^’^(R"*), let us set 
/3 = d{(f>dx^ A .. A da;"*) (this is a closed (m — l)-form) and a = 'ip. The explicit expression of da A is 
then dx^ A ... A da:"*. We can apply this to our situation for any couple (a, s). 
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endowed with the norm ||/||^i = sup \^t* f\ ■ The particular choice of 

0 <i<l Li(R4"i) 

does not affect the resulting space and gives an equivalent norm. 

The PDE (23) for a component of n on a geodesic ball Br C M can be rewritten 
in normal coordinates in the following form (implicitly summing over i,j and dropping 
the index [3 indicating the component for the rest of this section in order to keep a 
lighter notation) 


for f £ (set / = 0 in the complement of Br). Here (5®-^ denotes the Kronecker delta, 
A is the usual Laplacian and — 5®-^) = 0(|xp) and ^i{^/\^\g^^) = 0(|x|). Up 

to a dilation of the (small enough) geodesic ball Br C Ai we can assume that we are 
working on a ball Bji{0) with radius R much larger than 1 . We suitably re-define u, 
g and / replacing them by n 9 / (j^-)- By abuse of notation we will 

however keep the notation u, g and / on Bf(. With this in mind we have (-\/[^ 5 ®-^ — W) 
and ^i{^/\^\g^^) much smaller than 1 on Bii{0). 

Consider a smooth bump function x compactly supported in Bji{0) and identically 
equal to 1 on B^/ 2 {^)^ with |Vx| and |V^x| much smaller than 1. Extending u by 
setting it identically 0 in the complement of 5^(0), the function yu is well-defined on 
the whole of Compute 

^ixu) = Ax u + VxVn -|- = 

(Ax) n Vx • Vn - X 9 '^DijU - xt^Dju + X\/M/> 

where /r®-?' = and A = di{^/\^\g^^) are smooth functions of x G M^®®® and 

small in modulus. The latter PDE is valid on M^®®® thanks to the multiplication by x- It 
can be studied as a perturbation of the PDE A(xn) = X'v/l^/j thanks to the fact that 
the remaining terms are small. Recall that if v solves the PDE An = / with / G /t^ we 
have that v G and ll^^lliya.i < C||/||^ 1 . 

Consider the space F = {w £ g /t^(]R^®®®)} with the norm 

||tn||^i,i -|- ||V^t(;||^i. We will find by a fixed point technique a unique solution in 
F to the PDE 


I Av = {Ax)v + Vx-'^v-x9'^DijV-XT^DjV + X\A9\f (25) 

\ v{x) —>■ 0 as |x| —)• oo 

For this purpose, recall the Calderon-Zygmund operator CZ ■ f ^ f |j._^| 4 m f{y)dy- 
This operator is justified by writing the solution to An = / on M^®®® (the unique solu¬ 
tion under the condition “boundary data 0 at infinity”) as convolution with the Green 
function G{y) = '■ differentiating n twice yields that V^A“^/ is controlled by 

the convolution of / with the Calderon-Zygmund kernel on M^®®®, i.e. f f{y)dy- 

A deep result, see e.g. ED, shows that CZ is a bounded operator from to 

For any w £ F we can find a unique solution v £ F (with good estimates for the 
norm || • \\f) to the PDE 

Av = (Ax) n; -k Vx • Vw - X 9 "^DijW - xr^DjW + xVWlf, (26) 
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with boundary condition v{x) —)• 0 as |x| —)> oo. To see this, let us analyse the rdr.s. of 
(26). The function DijW belongs to since w £ F and multiplication with a smooth 
compactly supported function preserves the propertjj^ of being in . The same goes 
for X\/H/- The functions w and Vtc are respectively in and in W^’^, that are 
subspaces of Thus all the terms on the r.h.s. are in k^ and so is their sum. Writing 
equation (26) as Av = / we have seen that f £ k^ and the unique solution is given by 
f i^_j 4 m- 2 f(y)dy. In addition to ||u||^y 2 ,i < C\\f\\i,i we also have 


||V\|Ui<||C'Z(/)|Ui< Const ll/Li. 

Let us prove that the map w ^ Tw = v just constructed is a contraction from F into 
itself. Given wi,W 2 £ F, the difference vi — V 2 ■= Twi — Tw 2 solves 


A(ui -U2) = (Ax) {wi-W2)+Vx-'^{wi-W2)-xy-'^Dij{wi-W2)-XT^Dj{wi-W2) 

and by the estimates recalled above we can see that \\Twi —Tw 2 \\f is bounded by the 
Hardy-norm of the r.h.s. Recalling the smallness assumptions on the moduli of V^X) 
Vx, and we hnd that 


llTtci — Tw2\\f < e||rci — W2\\f 

for some e < 1, i.e. T is a contraction. The Banach-Caccioppoli hxed point theorem 
yields a unique hxed point of T in F, in other words a solution v £ F oi the PDE 
(25). We also have ||u||i7’ < e||u||i? -|- const||x\/H/||/;i, so absorbing the hrst term on 
the r.h.s. 


uIIf < Const||x\/H/||/;i- 


Equation (25) reduces on to 


A-c + - d'‘^)dijv + di{4/\g\g'-^)djV = \/H/, 


since x is identically 1 there. So v and u satisfy the same PDE on Bji/ 2 - Their difference 
thus solves 


di (^y/\g\g'^dj{u - u)) = 0 on Hr/ 2 . 

So by standard L^-theory for divergence type elliptic operators u — v is smooth and we 
have ||n — < ConstF\\u — Recall that u is to start with, 

since it was a coordinate function for a map taking values in the compact submanifold 
Af C M'^. Choosing p = > 1 (by the Sobolev embedding v £ which 

embeds in LP and ||u||iP < const||n||^ 2 ,i) and recalling that the LP^’^-norm is controlled 
from above by the T-norm, we have 

< ConstR {\\u\\LPi^Bii/A) + ^ ConstR {\\u\\lp{^Bh/a) + Il'*^llty 2 .i 


< ConstR 





+ 


xVWlf 


l£ 


< ConstR (||m||oo + ll/Li), 


^■^Let a £ F and h £ L°°. Then * {o-b)) (*) = / ^t{x — y)a{y)b{y)dy and so |$t * {ab)\{x) < 
||&||oo|'t>t * aK®), whence ||a6|hi < ||6||oo||a|hi • 
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from which it follows that 


< \\v\\f + \\u - < ConstR (||u||oo + ||/Li) • 

The geodesic balls Br in the beginning must be chosen so that - 6^/4 form a finite cover of 
Ai and — 5^^ is uniformly small in the C'^-norm (this is possible since At is compact 
and we have smooth data on it). Then since all the radii r are uniformly bounded away 
from 0 and oo the dilations that send Br to Br produce a fixed multiplicative constant 
when the final estimate is rewritten in the original ball Br- We can then add up these 
estimates as the geodesic ball ranges over the finite cover. 

Repeating argument above for every component of the map u : Ai ^ {M ^ 
we have the desired result 


u £ W‘^’^{Ai) and ||u||w 2 .i(a 4 ) < Const{^_M,g) + E{u)j , 

where E{u) denotes the Dirichlet energy of u. Remark that we have also shown that 
G with norm bounded by + E{u)^, therefore it also holds 

Au G (with A rather than A^) with the same estimate on ||Am||^i. 

3.3 Consequences: ^-regularity. 

The PDE Am G leads to Morrey decay estimates showing an e-regularity result 
for a triholomorphic map u G W^'‘^{Ai^J\r) satisfying (|^, i.e. we get the same results 
that are well-known in the case of stationary harmonic maps, although our triholomor¬ 
phic maps are only almost-stationary. The proof exploits the duality bmo — /t^, just 
as in the fundamental works on stationary harmonic maps |10) and |5]. Due to the 
importance of the result for the subsequent development of the paper and for the sake 
of clarity and completeness we give in this section the explicit argument for the key 
lemma in this regularity result. 

The space of “local bounded mean oscillation” is defined as 



(27) 

where fq denotes the average on a cube Q) /q = |^ / f{x)- The space bmo is the dual 

Q 

of see m 2.1 .2] and [38l 2.3.5] (|T2] for the homogeneous spaces and BMO) 
with duality pairing 


(/> 9) 


fg for / G bmo and g £ H}, 



ll/llbmoll^Li. 


We can then obtain Morrey’s Dirichlet growth lemma under a small energy assump¬ 
tion. 


Lemma 3.1. There exists eg > 0, r > 0 and rg > 0 depending only on the geometric 
data on At and Af such that whenever u £ W^’‘^{Ai,N) satisfies a.e. the triholomorphic 
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equation 0 and the local strong approximability condition and there exists Br{x) such 
that r < ro and |VmP < Sq, then it holds 


1 


(rr) 


4m—2 


Srr { 2 ^) 


|Vnp < - 


2 r 


1 

4m —2 



Proof. This argument is due to |10) . where S^-valued stationary harmonic maps were 
studied. The argument goes though in our situation as well, with little modifications, 
thanks to the similar PDE structure and to the almost monotonicity of the energy 
ratio. Regarding the latter, remark that (by the compactness of Ai) there exist i?o > 
0 and C > 0 such that, for any x € A4 and whenever ri < r 2 < Ro, we have 
f f |Vup. This will be needed in the final stages of 

Br^ (x) Br2 (x) 

the proof. 

We can assume that t < ^ since we need the result for some small r. We will argue 
by contradiction and thus consider a sequence 


Ua 


BrAXq) Af C with 

Tn 


1 


^4m —2 


|Vuqp = Ag —)• 0 and —)• 0 as g —)• 00 (A^ / 0 ) 


Brq i^q) 


and j / |Vu ,|2 for all g. 

Brrq{Xq) Brq{Xq) 

We begin by rescaling each map: denote by Oq := {uq)x^,rq '■= ^( 3,^)1 Jsr (xq) 

define 


Uq 


and 


Vq : Bi{0) 


S Vq{y) = 


Uq{Xq + rqy) - {Uq)^ 


Then 


/ iVuqP = 1 and / |VugP > - (28) 

JBx( 0) JBr(0) ^ 

for every q. Remark that the VqS have zero average, so they have a unform W^’^-bound. 
We can then extract a subsequence (not relabeled) such that 


Vq ^ V strongly in L^(Si(0)) and Vvq —>■ Vu weakly in L^(i?i(0)). 


The map v is then harmonic: this can be seen as follows. Consider the PDE (23) 
for Uq, substitute Uq{x) = AqVq{f- — Xq) + Uq and rewrite the PDE (changing variable 
y = — Xq) on i?i( 0 ) for Vq: the l.h.s. becomes 

r2 

q 

Let us see how the r.h.s. transforms. To avoid a heavy notation, remark that the 
terms that we must analyse are of the form d {p,{x)I{u)) du for smooth functions // and 
I (depending only on the geometric structures). Moreover each of these terms comes 
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with an associated term with which it forms a Jacobian. After the substitution we get 
for each such term 



and setting flq{y) = y{rqy + Xq) and multiplying by ^ the latter becomes 

{dyyq){y)I{XqVq{y))dyVq{y) + ilq{y){d^I){XqVq{y))dyVq{y)dyVq{y), 

and therefore the PDE for Vq has the schematic form 

AgVq = {Vflq)I{XqVq)Vvq + Xqflq {Vug, VVq}j + J-terms + AT-tenus, 

where {’Vvq,'Vvq}j denotes a Jacobian term and the other functions appearing are 
smooth, bounded and depend on g, i Vi, /, VI (similarly for the remaining J and 
AT-terms). Let us look at the hrst term on the r.h.s. Remark that V/ig —?■ 0 in 
as O' —)• oo (recall that Vq —?■ 0) whilst I{XqVq)Vvq are bounded in independently 
of q. Let us look at the second term on the r.h.s. Thanks to the Jacobian structure 
it is in the Hardy space /t^ with /t^-norm (and thus L^-norm as well) bounded by 
XqK iViJqp = XqK for a hxed constant AT. In the weak form the PDE reads, for 
every ip G C“(I?i(0)) 




ipO{rq)VVq + Xq 



ipjlq {Vvq,Vvq} j + J-terms + H-terms 


and so 




we conclude 


< (Aq+ rq)A'||(^||oo —t 0. By the weak convergence Vvq —)■ Vv 



{Vv,Vip) = 0 for all ip G C“(Hi(0)), 


i.e. V is harmonic on In particular by Bochner’s formula |Vt;p is subharmonic. 

Therefore the average of |Vup on Bt{0) is controlled from above by the average of 
|Vfp on l?i(0), i.e. 


[ \Vv\‘^<r^[ \Vv\^ = T^<l. (29) 

T ^ JBriO) Jbi{0) 4 


The proof will be concluded by showing that the convergence X/Vq 
strong in A^(I?i(0)), thereby reaching a contradiction with (28). 


—>■ Vf is actually 
This is the more 


delicate part of the proof and here the duality bmo — plays a fundamental role. 

Consider C, G C')T(I?i(0)) with 0 < C < 1 identically 1 on Hi(0) and identically 0 
on the complement of B5{0). The previous weak equations for Vq and v give, for all 
test functions ip, 



Vv,Vip) = Xq [ 

iBi(O) 


ip {{X7Vq,X/Vq}j + {X/Vq,X7Vq}j + {X7Vq,X7Vq }+ 

(30) 
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+ / 0{rq)ipVVq. 


We will now set ip = C'^{vk — v) in ordei^^ to obtain an L^-estimate on V{vk — v) as 
follows: the l.h.s. of ( |30| ) becomes 

/ Q^\V{Vk-v)\^ + 2 ! (^{y^-v){V{vk-v),VO> 

JBi{0) JBi{0) 

> [ | VK -^;)|2 + o ( 1 ), 

2 

by the strong convergence of fg — f to 0. It will thus suffice to prove that the right 
hand side of (30) goes to 0, in other words that 


/ - v) ({V?;g, Vr;g }j + {Vvq, Vvq}j + {Vvq, Vug}^) 

JBi{0) 

is bounded independently of q. We will regard this integral in the spirit of the bmo — fi^ 
duality. Recall the Hardy space control (independently of q) 


II {VVq,VVq]q ||/;l(_B^(0)) + || {Vl^g, Vfgjj |Ui(_B^(0)) + II {Vfg, VUg}^ Li(Bi(0)) < 

we will thus aim for a uniform ftmocontrol of C‘^{vq — v). For z G ^ 7 ( 0 ) and r < | we 
have by the Poincare and Holder inequalities 



The Morrey-Campanato norm that appeared on the r.h.s. is a natural quantity for 
stationary harmonic maps, it is invariant under dilations of the domain and we know 
that it has good monotonicity properties. Indeed, going back to Uq (from the definition 
of Vq) we have 


1 


r.4m—2 


I \^Vq\^ = 


1 


1 


( rgr )^™'“2 


I IVWgl" < 


Br{z) Brqr{Xq+rqZ) 

using the almost monotonicity, and then the inclusion Bi^ (xq + Vqz) C Bj.q{xq), 


1 8^™-2(l + Crg) 


< 




(fq) 


4m—2 


IVUgP < 


2 ^ 1 8“(1 + Crg) 


Bl (Xn+rqZ) 




(rg) 


4m—2 


|V«g |2 = 


Brq (^g) 


= 8“(1 + Crg) <C^,g. 

Recalling (31) we can see that Vq have equibounded 6mo-seminorms on i?7(0) (i.e. the 

8 j__| 

mean oscillation on small balls). In order to complete the control of the bmo-norm (27) 
it is enough to recall that the Vq^s are bounded in because they have zero average. 


^®This test function is in Wq’^{B i{0)) n L°°{Bi{0)) rather than in such test functions 

can be seen to be valid by a density argument. 
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Moreover by the John-Nirenberg inequality the Vq's are also equibounded in L'P{Bi (0)) 
for all 1 < p < oo (we will need later p = 4m). Consider now C,Vq^ we need to control, 
for z G and r < |, 

\bM\ / - \bM\ / 

Br{z) Br{z) Br{z) 

By the smoothness of C, we have {Cvq)z,r = Q{p){vq)z,r for some p G Br(z), and therefore 
the last term in the previous inequality is bounded by (here C denotes a constant 
depending only on C) Cr .^ Is, J \vq\. The first term on the r.h.s. on the other hand 

" ^ BrC) 

is controlled by J \vq — {vq)z^r\- Putting together and using Holder inequality 

" ^ Br(z) 


\Br{z)\ 


\CVq 


Br{z) 


{C'^q)z,r 


< CCm,g + 


c 

,^4m—1 




1 

4m 


^4m-l 


and these are equibounded in q by the L^^-estimate above. The estimate of the bmo- 
seminorm for 2 ; G is trivial since C vanishes on the complement of Bs^g^. We 

have thus found that the r.h.s. of (30) goes to 0, i.e. Vn„ —)• Vv strongly in L^(Bi (0), 
thereby providing a contradicion between (28) and (29). □ 


Iterating the previous result at smaller and smaller scales and following the standard 
De Giorgi/Morrey’s scheme one shows that whenever the energy density at x is below 
£0 the decay ^ 4^,2 f |VmP < Cr''' holds for some 0 < 7 < 1 , from which the uniform 

Br{x) 

^-Holder continuity of u in Br{x) follows (see |16) . |25) . |19)). The Holder continuity 
can then be bootstrapped to (7°°-regularity. This last step can be achieved by standard 
elliptic theory just as in the case of weakly harmonic maps. Indeed we will now show 
that, also in the case of almost harmonicity that we are dealing with, the PDE satisfies 
suitable assumptions. 

The continuity of u in B^(x) allows us to localize not only in the domain but also 
in the target, therefore we can write u : B^(x) —>■ Abusing local charts in Br{x) C Ai 
and around u{x) G Af (let us use coordinates respectively x^ and y^). Then taking a 
local variation of u in the target amounts to comparing the energy of u with that of 
u + t(j) for (j) = G C^iBr (x)). We can compute, as in Section (analogously 

to (Q with equality for t = 0 and the computations that followed) 


dt 


t=o 



\Viu + t(j))\^ = 


(32) 


= j d((a;0''"-')A(n*(r^H,)) + d((a;,f™-i)A(n*(.^Hj))+d((a;fc)2—i)A(n*(r^Hx)). 

Br{x) 

The l.h.s. of ( |32[ ) can be rewritten using the expression of the metrics g and h (respec¬ 
tively on Ai and Af) in the focal coordinates and leads, with standard computations 
(see e.g. |l9j or | 2 T|), to the expression 
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With the r.h.s. of (32) we are thus introducing a perturbation term in the weakly 


harmonic map equation. Let be expressed in local coordinates as a sum of terms of 
form Cf^u{y)dy^ A dy'^. Then u* is a sum of terms of the form Cfj,y{u)(j)^ A dv^^ 

therefore the perturbation term in the weak harmonic map equation is of the type 
f{x,u,Vu) with smooth dependence on the variables and linear dependence in the 
first derivatives of u: it satisfies therefore all the assumption (A1)-(G2) of |21[ Section 
6.2], The smoothness of u thus follows. 


We can further obtain the following result. When u is as in Lemma [3.1 [ there exist Sq 
and ro suchj^ that on each ball Br{x) with r < tq and on which ^ 4^,2 | Vnp < Eq 

we have that u is smooth on Br{x) and the following gradient estimate holds: 


sup |Vn(?/)| < 




y&Br{x) 


The proof of this estimate requires a blow up argument along the lines of 

2 . 2 ]. 

To conclude we have 


(33) 

Theorem 


Proposition 3. There exists Sq > 0 depending only on the geometric data on M and 
Af such that whenever u G satisfies ^ and ^ then 

0(ti, x) < Eq ^ X is a regular point of u, 

i.e. there exists r > 0 such that u is smooth in Br{x) (and actually Q{u,x) = Oj. 


Moreover (33) holds. 


(34) 


In view of this, the singular set of u, i.e. the closed set 

Sing„ := {x £ Ai : $r > 0 such that u is smooth in Br{x)}, 
can be characterized as follows 

Sing„ = {x £ M : &(u,x) > eo}- 
Proposition 4. With the same assumptions as in Proposition^ we have 
Sing.^ is a closed set and Bfi'^~'^{Singfij = 0. 

The proof is rather standard and involves a covering argument and the monotonicity 
formula, see e.g. |191 3.5]. (in our case we use the almost monotonicity instead but no 
change is needed). 


°We can of course take the £0 and ro here to be equal to those from Lemma 3.1 
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4 Blow up analysis 

Consider the following compactness issue. Let be a sequence of tri- 

holomorphic maps from {A4,g) into Af C with d('u|a) = 0 whenever a is a closed 
2 -form on Af. Assume a uniform energy bound, i.e. there exists A > 0 such that 
< A for all i. Then, up to the extraction of a subsequence, that we do not 
relabel, we have ug ^ u weakly in for some u G A/^. What can we say 

about the limiting map u! 

Recall the almost monotonicity result. Consider the set (here C , Eq and tq are as 
in Proposition Lemma |3.1| and Proposition recall that they depend only on the 
geometric data on AA and AP) 

T, := < X G Af such that 

0<r<ro ^ 


lim inf 

i^CSD 


1 -|- (4m — 2)r 


r>4m—2 


(a^-^) A u*ni + (ar“ ) A u*nj + A u*nK 


> £0 


= n 

0<r<ro 


< X G A4 : lim inf 

^oo 



Br{x) 


|Vtt£p > £0 


Cr . 


(35) 


This is called the blow-up set, or bubbling set: it is the (closed) set where the energy 
concentrates and 

Fact: away from S we have locally smooth convergence of U£ to u. 

This can be seen as follows: consider x ^ S, i.e. there exists 0 < r < ro such 
that for a subsequence ui- we have f < Sq and thus u^. are smooth 

Br{x) 


in Br{x) and satisfy (33) in i?r(x), so that we have ’“-convergence of U£. to u in 
Br(^x). In particular u is C^’" on AA \S. Remark that by the (7^’“-convergence we can 
pass the PDE to the limit and ?x is a triholomorphic smooth map on Br{x). The local 
C^’“-convergence can then be bootstrapped to C^-convergence for any k. Moreover 
^ on A4 \S as Radon measures. Following a well-known 

scheme involving a covering argument and the almost monotonicity of the energy ratio 
(see e.g. |23l Lemma 1.5]) we have 

Fact: S has locally hnite ^-measure (and locally hnite Minkowski content). 


This is proved by the following local argument. Take SnRi, it is compact. Consider 
0 < h < and take a hnite cover of SCRi with balls Bs{xj) of radius 6 and centered 
at points xj G S. By Vitali’s covering result we can choose a subfamily of these balls 
such that B^{xj) are disjoint and B'^s{xj) covers S H Ri. For all ^ large enough we 
have that (for every j) IbsIx^) -C6 > ^. Therefore 


54 m -2 < 

j 


£0 JujBsixj) 


|V«,|2 < - / |Vn^|2 < -A, 
£0 Jm ^0 


where we used that the BsixjYs are disjoint. Then we have, recalling S n Ri C 
UjB^sixj), that 
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^0 

Remark that this result can be refined, in the following way. Since UjB^sixj) covers 

S n Bi, then the open set UjB^sixj) contains the 5-neighbourhood of S n Ri, i.e. the 

set {y : Bi : dist(y, SnRi) < 5}, therefore this 5-neighbourhood has (4m)-dimensional 

measure at most 4^™'^A5^. In other words the 2-dimensional Minkowski content of 
_ eo 

S n Ri is bounded by 4^™^ A. 

By ’“-convergence to tt on A4 \S, which implies the strong PR ^’^-convergence, 
we also have, for any 2 -form a on A/", the strong L^-convergence tt|(a) —)• u*{a) when 
we restrict to an open set in the complement of S. Therefore u satisfies the strong 
approximability condition away from S. However this does not necessarily hold across 
S, where the PP^^’^-convergence will generally fail to be strong. We now analyse this 
aspect further. 

Consider the Radon measures | Without loss of generality we may as¬ 
sume \'VUil'^d'H‘^^ /4 weakly as Radon measures. By Fatou’s lemma y = | 

V for some non-negative Radon measure i/ on AP. The previous discussion then gives 
that v is supported inside S and Sing„ is a subset of S. On the other hand, whenever 
X € Yl then we have that for every 0 < r < tq there exists L G N such that and for 
every i > L it holds f |Vn£p > ^ and therefore for a.e. 0 < r < rg we must 

Br(x) 

have If w is smooth at x then ^ 4^,2 / |VmP = O(r^) and therefore for 

Br(x) 

a.e. small r we have i.e. x is in the support of i'. We conclude therefore 

Fact: S = Sing„ U spti^. 


Similarly to |23) Lemma 1.6 we further have the following facts: 

(a) For every x (z M the ratio — im -2 almost monotone in r, i.e. for r < rg 

it is the sum /(r) -|- 0(rf(r)) of a non-decreasing function of r and an infinitesimal 
function. This allows to define the density of /r for all x, 0(/r, x) := limr-s-o ■ 

This follows by passing to the limit — >■ yL{Br{x)) for a.e. r < xq and by 

using the almost monotonicity formula for U£. 

(b) X G S 0(^,x) > Eq. Indeed, if 0(/i,x) > eg then for r < rg we have 

> £0 -Cr by (a) and therefore liminf^^oo ^ 3^=7 |Vu£p > eg -Cr by the 

weak convergence of measures and thus x G S. On the other hand if x G S then for 
a.e. r < rg we have > eg —Cr and sending r —?• 0 we find that 0 (/r, x) > eg. 

(c) For the PP^^’^ function u the set |x : limsupp_,.g ^ 4^,2 |Vnp > o| = {x : 

&u(x) > 0} is a set of zero TP^^^^-measure, see |11| 2.4.3. and 4.8]. In the complement 
of it u is approximately continuous (in the sense of Lebesgue points). Moreover we can 
pass to the limit, away from E, the second order PDF that we obtained in (32) because 
we have (7^’“-convergence, i.e. on the complement of S the limiting map u satisfies the 
weak formulation of 


AgU + A{'Vu,'Vu) = f{x,u,'Vu), ( 36 ) 

that is (in a first moment we must take cf) G C^{Ai \E)) 


25 













J VuV(p + A{Vu,Vu)ip = j f{x,u,Vu)ip on Ai, 


(37) 


where A denotes the second fnndamental form of A7 C and / is a smooth fnnction 
depending only on the geometric strnctnre on M and J\f and depending linearly on the 
third variable. Moreover this PDE extends across the 7 ^^”^“^-dimensional set S by a 
standard capacity argnment, so u satishes the weak formnlation of this second order 
PDE globally, i.e. we can allow (p G in (37). At a continnity point of u we 


can localize it in the target and nse the resnlts from 
smoothness of u. 


as we did in (32) to infer the 


/4m—2 


G S the map u satishes Q{u,x) = 0 and it is there 


Fact: for "-a.e. 

approximately continnons. 

The three facts (a)-(b)-(c) imply that the density of /i at points of S and the density 
of v w.r.t. the measnre coincide a.e. on S (in particnlar the density of u exists) 

and the valne of the density is > 69 . From [271 we dednce 

Fact: the set S is (4m — 2)-rectihabl^^ The measnre v can be expressed as 
u = 0 (x)LS and 0 (x) > Eq for a.e. x. 

The limiting map u solves and (37), however the weak IT^’^-convergence Ui ^ u 


is not enongh to infer that d{u*a) = 0 whenever a is a closed 2-form on N. This canses 
a possible lack of almost-stationarity of tt, which is a very important and possibly hard 
issne to analys^^ Focnsing on the case in which A4 is hyperKahler (i.e. the maps ui 
are stationary harmonic), all we can say abont the limit is that the conple (u,Ti) is 
stationary in the sense of |22|P^ bnt a priori neither S nor u need to be. Determining 
whether they are is an open qnestion, as well as the regnlarity properties of both S 
and u. 


5 Quantization of the energy - proof of Theorem 1.2 


1.2 


Consider U£ : Ai ^ Af, a seqnence of W 


1,2 


In this section we prove Theorem 

triholomorphic maps satisfying (|^ with a nniform energy bonnd f \Vui\‘^ < A < 00 for 
all ^ G N. By the resnlts in the previons section we have U£ ^ u weakly in where 

u is a triholomorphic map, with strong convergence away from a set S that is (4m — 2)- 
rectihable and has locally hnite 'H^™'~^-measnre. The map u is smooth away from the 
set Sing„ C S, and is approximately continnons away from a P^^™’“^-neghgeable set. 
We moreover have 


with 0 > eo > 0 a.e. on S. The qnestion that we analyse in this section is the following. 
Is it trne that for x G S we have 

e{x) = f^E{cp,), 

S = 1 

|23j Lin gives a self-contained argument for the rectifiability, without relying on the deep 
recifiability result of Preiss. 

^®This was overlooked in j^. 

^®Roughly speaking, the first variations of E as a varifold and of u as a map might fail to vanish 
but they annihilate each other. 
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where each : 5^ —)• AA is a (smooth) non-constant harmonic map? This is known 
as quantization of the defect measure. We will answer the question affirmatively. In 
the case that M and M are hyperKahler the affirmative answer was already given in 
m with a different approach. We will give an alternative proof, covering also the the 
more general case we are analysing, with Ai almost hyper-Hermitian. Our proof relies 
on Lorentz space estimates, as in |24| . exploiting the structure of the PDE that we 
analysed in Proposition and the W^’^-regularity result from Theorem 0 In view 
of Section we will go through the arguments quite comprehensively, omitting some 
details for which we refer to |24) . 


For the purposes of Theorem 1.2 it is enough to restrict our attention to points 
X £ T, where the approximate tangent TxT, exists and such that the limit map u 
is approximately continuous at x: by doing so we only neglect a P^^^^^-negligeable 
subset of S, which does not affect the result. With this in mind, at any such x and for 
every i we can find a suitable dilation so that (using a diagonal argument) the 
dilated maps, that we will denote by U£, converge weakly in to the constant u(x) 
(the value taken at x by the precise representative of u) with blow up set T 2 ;S|^ We 
thus have ^ and the measure iz is a multiple (with factor Q(x)) of the 

^-dimensional Haussdorf measure restricted to the linear (4m — 2)-dimensional 
subspace 


In view of the quantization question that we are addressing, we can therefore work 
straight away on the U£ and prove that 0(x) (that is now a constant on is a sum 

of energies of non-constant harmonic spheres. For the sake of notational convenience, 
we will drop the tilde for the rest of the section and assume that is x (0,0) 

in suitable coordinates. So we assume that we are dealing with a sequence of strongly 
approximable triholomorphic maps U£ : i?^™'“^(0) x Bf{0) —?• A? with uniformly 
bounded energies such that 


U£ cnst weakly in ^(0) x 

where cnst stands for a constant map and the convergence is strong (even C^’") locally 
in the complement of i?^™'“^(0) x (0,0); moreover 

\Vu£\'^d'H^^ 0 L {B^^-^{0) x (0, 0)) for some 0 < Eq < © e K- (38) 


From the monotonicity applied to U£ at the points 0 and Ci, ^ 2 ) •••; ? 4 m -2 (where the 
latter (4m — 2)-points span x (0,0)) we have, as in [23l (2.11)], that 


4m—2 

E 


j4m-2(o)x52(o) 


a=l 


duf 


dXa 


0 as £ —oo. 


With the notation Xi = (xi, ...,^ 4 ^- 2 ) and X 2 = ix 4 m-i,X 4 rn) we then define 


/ 4m—2 

E 


du£ 


dxn 


2 

iXi,X2)dX2 


^°This is a special case of Federer’s reduction argument, that is performed at all points in I23h for 
a general exposition of this argument we refer to |33| . Appendix A. 
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and by Fubini’s theorem we then have 0 in ^(0))- Recall the Hardy- 

Littlewood maximal fnnction Mf^ of fe 


Mfe{Y):= sup 

0<r<i ' 


fe{Xi)dXi for y G bT-\o) 


Br 


Hy) 


and the weak-L^ estimate 


n 


4m—2 


{{Mfe>X})<^\\fi\\Li 


So for every (i > 0 we can find C Bf^-^{0) with {Bf^-'^(0)) - 

2 2 

Cmd and 

4 m -2 / ^ for ^11 Y ^ As and for all large enough £. 

r" J 


0<r< 


Br 


XY) 


set 


Recall further that for each £ the map Ui is smooth away from a T^^^^^-negligeable 
Taking a countable union we conclude that for 'H^”^“^-a.e. Xi G B‘^^~‘^(0) and 
for all X 2 G R^(0) all maps Ui are smooth at {Xi,X 2 ). Therefore we have a sequence 
of points Xl G such that 


each map U£ is smooth near {X\} x Ri(0) 
and sup j /f(Xi)dXi0 as ^cx). 


(39) 




(Xi) 


Now we need to “track down the first bubble”. For this purpose we show that, for 
every £ large enough, we may find 5^ G (O, and G B\ (0) such that 


max 

X2&Bl^^(0) 5\ 


4m—2 
£ 


\Xui\'^dx 


(40) 


'(X()xB| (X 2 ) 


is achieved at X 2 with value -and X 2 —?■ (0, 0) for £ —)■ 00 . 

8 ‘ 2 


Indeed, the smoothness of - first line of (39)- gives that for each £ there is a 6{£) > 0 
such that for all X 2 G B\ (0) we have 


1 


lj4m —2 


\Xui\-^dx < - 


1 £0 


2 8 • 2-^^ ’ 


(41) 


Bl”--Hxi)xBj{X2) 

whenever 5 < d{£). On the other hand we may fix any arbitrary 5 > 0 and have, for 
all £ large enough, that 


max —-^ 

X2eBi^^(o) ( 54 ”" 2 


\Vue\‘^dx > eg. 


(42) 


j^Am-2(^l 


(Xf)xi?|(X2) 


are using Propositionand the almost continuous representative of a map as we did in 

section 
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for if that were not the case, we would use the e-regularity result Proposition 
and conclude that \Vu(\ < on x This would give (7^’“- 

convergence of Ui on x - 6 ^ 4 ( 0 ), contradicting the existence of the blow up 


set “^(0) X (0, 0). Conditions (41) and (42) imply (with the smoothness of which 
gives continuity in 5 of the energy ratio at a point) that for every large enough i we can 
hnd 61 satisfying (40) with —>■ 0 as £ — 00 . Finally let us check that X 2 —t (0, 0). If 


that were not the case, i.e. if there exists c > 0 s.t. I^H > c for all £ in a subsequence 
(not relabeled), then we would have (using the almost monotonicity of the energy ratio 
for balls centered at {X^^X^)) for all I large enough so that c > 45^ 


\Vu(\^dx > Cc ,m,£Q 




'(Xf)x(B2 (x|)\B2 (X|)) 


independently of i for all large enough 1. However this contradicts the strong conver¬ 
gence un —)> cast on ^^’"“^(X^) x (^ 4 ^ 2 (-^ 2 ) \ ^ 012 (^ 2 ))- 

We consider now the maps V£ dehned by V£{y) = U£ (^(Xf,X 2 ) + d£y). The map V£ 


is dehned on ^ x Bi^ := B"^^ ^ x and R£ —>■ 00 . Then we have 


V£ 


V in X m2) as i 


00 , 


(43) 


where v{y) = v{y 4 rn-i,yim) is a smooth nonconstant harmonic map from M^ into J\f 
with hnite energy (n does not depend on the remaining (4m — 2) variables in view of 
the second line of (39)). For the proof of the Cj^’^-convergence we refer to |231 page 


818]. By Sacks-Uhlenbeck’s result |32| we can see v as a smooth harmonic map from 
into Af (using the conformal change of coordinates induced by stereographic projection 
from the north pole). This gives the hrst bubble 

Remark that the ’“-convergence is only local on M^, therefore after the conformal 
change of variable that lets us see the convergence on we have a sequence of domains 
exhausting and leaving out a small disk (shrinking as £ ^ 00 ) around the north pole: 
the -convergence holds on any domain that leaves out a hxed disk around the north 
pole. On the remaining neck we need to perform further analysis. First of all we must 
locate other points where the energy is concentrating (locate other bubble domains) 
and where we can produce more bubbles by repeating the procedure illustrated for the 
hrst bubble: “dilation of disks when needed/reparametrization on 52” . We will hnd 
possibly several but surely hnitely many bubbles, since each bubble contributes at least 
with a hxed positive amount to the energy and we have a global bound. The second 
task will be to check if there is energy that accumulates in the necks connecting the 
bubble domains: this will be ruled out and will lead to Theorem Ol 


By Theorem 
now that 


1.1 


we have that ||Vn£||^yi,i < independently of i. Recall 

We recall the 


Pjembeds continuously in the Lorentz space L‘^dtTa >2 
notion of Lorentz spaces 


L 


2 , 1 / 


^) = 


L 


2,00 ( 


/ : II/IIl2.i := ^ VC"^{{\f\ > t})dt < ooj . 

{/ : ||/|li 2 ,oo := supt2/:”2({|/| > t}) < ooj . 

I i>0 J 
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Therefore, by first choosing good slices and then embedding, we get that there exist a 

constant A* and a snbset Ee C with measnre 

2 2 

snch that for all Xi G Ei we have 


||(Vn,)(Xi,.)||i2.i(B|(o))<A*. 


This choice of slices, i.e. valnes Xi G -B^, shonld be combined with the choice of 
El C with measnre C‘^^~‘^{E() > snch that the two conditions 


in (39) are satisfied for all Xi G E^. For the seqnel we will work on “good slices”, i.e. 


good valnes of Xi, meaning that Xi G H E^. We therefore consider the maps 
Ui{X 2 ) '■= Ui{X[, X 2 ) for a seqnence X[ snch thall^for every i we have 


each map is smooth near {Xf} x Bf{0) 

1 


0 <r<i 


j fiiXi)dXi ^0 asi^oo. 


||(Vu,)(Xf,.)||LM(B|(o))<AA 

From the blow up analysis we know that 


(44) 


0 = lim 

£^00 


|Vn,|2 


^V 2 ( 0 ) 


and the content of Theorem 1.2 is to show that this is a snm X)f=i ^{4's), where each 
4>s : —)• A/" is a (smooth) harmonic map. We have seen before, nsing (40) to find 


snitable points X 2 and snitable dilating factors 6 i, that V£(X 2 ) = Ui{X 2 + di{X 2 — X 2 )) 
converge strongly in to v harmonic from into J\f, which becomes the first bnbble 
(j)i after composing with stereographic projection from the north pole. 

For any fixed R the strong local convergence on Bji{0) gives | —?• | Vnp 

and by conformal invariance nix^) Ibr So the energy left (after 

the first bnbble has been acconnted for) is 


lim lim / 

R^oo 1^00 J Bl^^( 0 )\Bs^R(Xi) 


(45) 


(here V stands for DX 2 ) ^-nd we want to prove that it eqnals 'Ylis =2 ^{4's)- 
Observe that, for any fixed c > 1 and for any i? > 0 we have /g r\Br 

by conformal invariance ^ Ib^r\Br 

as £ —>■ 00 . Therefore 


lim lim / |Vtt£|^ = lim / |Vn|^ = 0. 

£^00 JB,s,R{Xi)\Bs,R(Xi) Jb,r\Br 


view of Section it should be stressed that we have a huge freedom in choosing this se¬ 
quence X}, namely for every ^ we only need to avoid a 'H‘*"*~^-small set, with measure at most 1% of 
£4m-2j'g4m-2^^ Sectionl^the 2-planes — X\ will be read in the original context (prior to the 


dimension reduction, i.e. prior to the transition from ui to ud). 
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Moreover we have for any c > 1 (remark that for all i large enough the domain 
-Bi (X|) \ Bj_(X 2 ) stays well away from S, since —)• (0,0)): 

2 2c 


lim / |Vn^|^ = 0. 

(X|) 

2 2c 


23 


Consider now the change of variables (conformal in X 2 G M^) induced by 

= U£{Xi,e~^,6) with X 2 = (r, 0) in polar coordinates. (46) 
Then j ’^2 {q)\b^ (x^) coincides with 

E{We,Be):= j iViT^I^ 

{Xf}x[log2|,log5iiJ|]x§i 

(where V stands for the standard gradient in t and 6 ) and we want to show that this 
is a sum of energies of harmonic spheres arising in the blow-up analysis; we know that 
this energy cannot accumulate close to the boundaries of the domain Bi, namely for 
any M > 0 (here c = with c used above) 

E{Wi,Q{0,M))-.= [ IVIT^P^O as^^oo 

J [log 2,log 2+M] xS^ 

lim Wu). E{Wi,Q{l,M)) ■.= lim lim / \VWe\^ 

R -^00 e^Qo R -^00 i^oo ^[1 iog5^_R|-M,| logi^RIlxSi 

Recall that there is a finite number of bubbles only. At this point we face a di¬ 
chotomy: either there is only one bubble or we can locate another bubble domain 
(in the new variables). If the second option happens, then we select the next bubble 
domain by choosing, for each a point G [log2, | \og5(^R\ — 1] such that 

Observe that 

lim [ |Vu^|^ = lim j = 0, 

{x|)\Bjl{x|) 1 (0)\S 1 (0) 

2 2c 25J 2cSJ 

i.e. the energies of vt, that actually do escape partly to infinity when extra bubbles are produced, 
must escape “in a controlled way”, i.e. the domains where these energies are localized (further bubble 
domains) escape to infinity slower than the annuli Bj_(0) \ B 1 (0) (for any choice of c > 1), i.e. 

25 ^ 2cS^ 

with a speed that is of lower-order compared to the speed of growth of the domains B 1 (0). We 

stress, in order to avoid confusion for the reader that is used to Parker’s bubble tree construction (26], 
that here (following [23], [24]) the bubble map produced by the conformal wrapping of the 2-disks 
(domains of the vi) onto is constructed, although equivalently in the end, differently than in |26| : 
here we are covering except for a small disk around the north pole (this missing disks shrinks to 
a point as £ —>■ 00 ) and the extra bubbling domains are escaping towards the north pole as well (at a 
slower speed); in |26| the bubbling domains are confined in the southern hempisphere and it is on this 
hemisphere that further bubbling analysis is performed. The same goes for the neck domain: here it is 
an annular domain with one boundary circle coinciding with the boundary of the missing disk around 
the north pole (thus going to the north pole as £ —>■ 00 ), while the other boundary circle is fixed when 
R is fixed (it comes from dBj{(0)) and goes to the north pole as —>■ 00 . In |26| . on the other hand, 
the neck domain is located close to the south pole. The difference originates from a different choice of 
the dilating factors S£. 
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^{Am—2)t 


rt+l 


|VW,|' 


|Xl|<e- 


Si 


(here V = (e ^Dxi, § 1 , -§ 0 )) achieves its maximum (for t G [log2, | log(5^/21 — 1]) at 


t = ti and this maximum is > £1 > 0 for some^'^ fixed £ 1 . Note that the map is 
triholomorphic and almost stationary harmonic with respect to V = {e~^Dxi, §ii^)) 
compare |24l (2. 12)0 Keeping in mind that the cyclinder [log 2, | logci^i?!] x gets 
longer and longer as £ —)• 00 , we find from estimates (47) that (respectively) ti —?• 
00 and I logfifi?! — —?• 00 : in other words, we can center the cyclinders at t£ and 
we will have that this translated cylinders get longer with lengths going to 00 on 
both sides (with respect to their centre) as £ —t' 00 . The centering of the cyclinders 
corresponds to the change of variable t = t — With the extra change of variable 
Xi = e~^^Yi we define 6) = We{Xi,t — t£, 9). We then have that Vi is defined 


’^”^■2(0) X [-Ml, Mi] X §1 with Ml 


on B. 

stationary harmonic (compare 


00 and Vi is triholomorphic and almost 


(3.3)]) and 


ini<i 



|VK,|2>£i, 


(48) 


si 


where V = The Dirichlet energies of Vi are uniformly bounded and 

we may assume that Vi —?• Voo weakly in (possibly strongly). The situation now 

is rather similar to the one that we were facing with the maps Vi when we produced 
the first bubble. The map Koo is smooth and harmonic on x (— 00 , 00 ) x 

and independent of the first (4m — 2) variables. We have two cases (compare bottom 
pictures of Figure]^ in Section]^: 

(i) If the convergence holds strongly in then the map Koo is not constant 

thanks to ( |48[ ). As such, by conformally changing its domain (— 00 , 00 ) x to S'^ \ 
{the two poles} we produce the second bubble (p 2 ■ S'^ ^ J\f (this is a case where further 
dilation is not needed and we can reparametrise on 5^ straight away). Intuitively we 
have replaced the piece [0,1] x with a sphere without two small disks. So instead 
of the original unique cyclinder (the neck domain we started from) we have a spherical 
domain (without two disks) where convergence of the Vi happens strongly, actually 

and two neck domains left from breaking in two the original cyclinder. On each 
of these two neck domains we start over with the bubbling analysis. 

(ii) If the convergence is weak then we proceed as in the case of (pi: we find a 

(non-empty, by (48)) blow-up set S, that is of the form x (t,9) and proceed to 

suitable dilations around concentration point of the energy (for a sequence of points 
converging to (t, 6)). We thus replace a small cap on the cyclinder (around (t, 6), this is 
the new bubble domain) with a new 5^ \ {a small cap} and attach them by a cyclinder 
(a new neck domain). The spherical part becomes the domain of <1)2■ On the new 
neck domain we start over with the bubbling analysis. Moreover, on the remaining 
part of the previous cyclinder (from which we removed a small disk) we proceed to the 
identification of further bubbling domains, if any are present. 


we can make sure that for any £1 > 0 this energy is below £1 then we must be in the case where 
we have only the Hrst bubble <j)i, as will be clear from the subsequent argument. 

^®There is a little misprint, the right hand side of [241 2.12] should read —e~^^AxiWi, where Axi 
is the Laplacian in the Xi coordinates. 
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This procedure will stop at some point because only a finite number of bubbles 
can appear, since each bubble contributes with at least a fixed amount of energy and 


this cannot become infinite (0 is finite). Recalling (45), now that all the 4>s have been 
identified, the next task is to show that the Dirichlet energies on the (finitely many) 
neck domains go to 0 as oo. To fix strategy and terminology (since topologically 
the necks are always annuli no matter how many times we perform a conformal change 
into a cylinder) we will look for an L^’°°-estimate for in each neck domain (thinking 
of it as an annulus on a sphere, as in the original picture obtained when we produced 
(j)i) by looking for an L°°-estimate in the elongated neck, i.e. the long cylinder obtained 
by the conformal change of variable ( [4^ . We conclude the section by explaining why 
we look for such an (apparently) different estimate and how we can get it. 

The Dirichlet energy is conformally invariant, so it does not matter how many times 
we conformally change a neck, therefore we can assume that the elongated neck we are 
looking at admits no further bubble domains. The aim is to conlcude that \Vui\ goes 
to 0 in the L^’°°-sense in each neck, then we will conclude that the L^-norm must go 
to 0 as well by interpolation with the uniform L^’^-bound. The L^’°°-estimate follows 
from a control of the type 


\X 2 — X^WX on each neck domain 


(49) 


for e as small as we like as long as we choose ^ and R large enough. Indeed (49) implies 


that whenever t > 0 and \Vx 2 Ue{X 2 )\ > t then we must have \X 2 — X!^ < This 

forces the set {X 2 : \Vx 2 U£iX 2 )\ > t} to be contained in the ball of radius ^ centered 
at X 2 . In particular the L^’“-norm of iVu^l is below e by definition. Now recall from 


(46) that 


r = e *, — logr = t. 


—dr d d 


and so \X 2 — X 2 \XX 2 '^ is nothing but the derivative of Wi after the conformal change 


of variables. In other words (49) amounts to a pointwise estimate on XWi in the 


elongated neck. Remark now that the conformal change of variables makes Wi almost 
stationary harmonic, so in every fixed size domain contained in a neck domain we have 
that the energy is as small as we wish in the L? sense (otherwise there would be another 
bubble domain there, as explained earlier) and by the interior gradient estimate we can 
conclude a pointwise L°°-estimate (as small as we wish) on VIT^. Once we achieve that 
iVu^l goes to 0 in the L^’°°-sense in the neck domains, using the uniform estimate on 
the L^d_j;^orms from (44) we conclude by interpolationp^ that the Lp' norm in the neck 
domains goes to 0, i.e. in (45) we only find energies of bubbles. We refer to m for the 
technical adjustments of these ideas. 

We conlude the section with the following 


Proposition 5. 

exists {ax,bx,Cx) with + bf 


With the assumptions of Theorem 1.3, for ^ 


-I- 

X ' 


-a.e. X G S there 
^ = 1 s.t. all the bubbles (fs ■ ^ M at x (using the 

notation in (^p ) induce holomorphic S‘^ ’s in M with respect to the complex structure 
— {uxl + bxJ + CxK). Moreover the approximate tangent plane at x toT, is holomorphic 
with respect to + bxj + Cxk. 


r2.1 


^One way to see this is by exploiting the duality, indeed for any v such that Vw € 

C M C we have /|Vnp ^ f Vv ■ Vv < || Vn|R 2 ,i || V«|T 2 ,=o. 


33 

















This is shown in m Lemma 2.2; we want to make transparent the important 
underlying idea and stress an aspect that will be useful to keep in mind in view of the 
arguments in Section]^ Let a: be a “good point” on S, i.e. exists and the limiting 
map u is continuous at x. The creation of a bubble at x is essentially due to the fact 
that the Dirichlet energy is accumulating at x: more precisely, the components VX 2 W 
are producing this concentration (here denotes the two directions normal to T^Tj). 
A suitable reparametrization in the normal directions X 2 provides the bubbles, whilst 
in the directions tangential to S the components Xx^ue. are going to 0 in This 
causes the bubbles to depend only on the normal coordinates X 2 . The equation Q 
passes therefore to the limit to each bubble with a degeneration in the Xi-directions, 
i.e. only the two derivatives in the directions X 2 are non-zero. It is easy to check that 
(i) any two directions identify a complex 2-dimensional plane with respect to a complex 
structure Qxi + bxj + Cxk for a certain choice of (ox, hx, Cx) G 5"^ (ii) the (4m — 2)-plane 
TxS is holomorphic for the same structure + hxj + Cxk that makes the normal 
2-plane holomorphic. A short computation shows that whenever we annihilate two 
directions in Q we find that the map becomes holomorphic with respect to ((a^i + 
bxj + Cxk), —{{axl + bxJ + CxK)), hence the result. The normal 2-plane can then be 
conformally transformed into an = CP^, as in the standard definition of bubble. 

Remark 5.1. It is worthwile noting the following: it cannot happen that at a fixed 
point X we produce some bubbles that are holomorphic with respect to (a^i + bxj + 
Cxk),—{{axl + hxJ + CxK)) and some others that are holomorphic with respect to 
{ttxi + bxj + Cxk), {{axl + bxJ + CxK)). With our choice of the PDE Q we are allowing, 
say, {i, —/)-holomorphic bubbles to appear but it can be checked immediately that if 
a map is triholomorphic as in Q and it is in addition {i, /)-holomorphic, then it must 
have rank 4 at least, whilst a bubble has rank 2. This feature is not typical of the 
general bubbling issue for stationary harmonic maps, where it could happen indeed 
that at a certain point we have some bubbles that are holomorphic and some others 
that are anti-holomorphic. 


6 Holomorphicity properties of blow-up sets - proof of 
Theorem 11.31 


Under the asumption that u is smooth on we know that S coincides with the 
support of the defect measure u in and we are assuming that this support is non¬ 

empty. We will prove the statement of Theorem 1.3 when £ is a Lipschitzian graph 
over some chosen (4m — 2)-plane. Let us see how the proof of the general case follows 
from this. For £ as in the assumptions of Theorem 1.3 we observe that for every point 
on £ we have a neighbourhood in which £ is a Lipschitzian graph. Consider the set 


£* = |x G S n : 3 an open neighbourhood Ax of x and ai + bj + ck with (a, b, c) G 5^ 

s.t. S n Ax = £ n Ax and this is smooth and pseudoholomorphic w.r.t. ai + bj -|- ck} . 

The set B is open in £ (also observe that the points in B are of density 1 for S). Let 
us show that B is closed in £. Since S is closed, we have that the closure in £ of B, 
denoted B, is contained in S. If there exists y £ B \ B we choose an open ball B' 
centered at y such that £ n i?' is a Lipschitz graph and remark that TiDB' has strictly 
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positive measure (it contains points of density 1). Then applying the result again, we 
obtain that in a neighbourhood of y the set S coincides with C and it is a smooth 
pseudoholomorphic submanifold, hence y £ V. Having now established that is 

open and closed in C, by the connectedness of £ we get that S n = C Cl and it 
is a smooth submanifold. Moreover locally it is pseudo holomorphic for a hxed almost 
complex structure, which gives a well dehned map from £ into S‘jj^ and this map is 
locally constant: the connectedness of £ yields that it is globally constant, completing 
the proof. 


In the rest of the section we prove Theorem 1.3 under the further assumption that 
£ is a Lipschitizian graph over a (4m — 2) dimensional plane. Recall that by the 


?4m 


rectihability of S proved in Section we have that for 'H'^^ ^-a.e. x G S H B^"‘ we 
have that x has density 1 in S and that there exists an approximate tangent that 


agrees with the tangent to £ at x. Moreover £\E is open in £. We obtain Theorem 1.2 


for all these (almost all) points on E H when we used the dimension reduction in 

the beginning of Section we had to ensure that at the point under consideration the 
limit map u was continuous in average and the blow-up set possessed an approximate 
tangent; the second condition is met now everywhere on E H We thus have, 

from Proposition!^ a well-dehned map from a.e. E n into the S'^ representing the 
complex structures al+bJ+cK on the hyperKahler M: it associates to a.e. x £ Eni?^™’ 
the complex structure for which the bubbles at x are holomorphic. 

Consider two points P and Q on ECR^”* where the approximate tangent to E exists 
and agrees with the tangent to £, obtained by the differentiability of £ (so at these 
points, that only exclude a null set, we have the quantization result). The strategy 
is to show that the sum of the bubbles at P and the sum of the bubbles at Q are in 
the same homology class in H 2 {N','Z). When this will be achieved, we will conlcude 
by means of a calibration argument. Recall that on M (with its metric h) we have an 
S'^-family of calibrations, i.e. the Kahler forms Ci(a,b,c) ■= + cQ.k for any 

choice of (a, b, c) with a'^ pb"^ = 1 and a 2-dimensional integral current is said to be 
calibrated by a(^a,b,c) when a.e. approximate oriented tangent plane T is such that the 
calibration <yi^a,b,c) and the natural area element induced from h have the same action 
on T. For the Kahler forms the notion of being calibrated by a^a,b,c) is equivalent to 
being holomorphic with respect to the associated complex structure al -|- 6 J -|- cK (see 
e.g. O page 4]). 

We need to recall at this stage, from Section how the analysis of the bubbling at 
P was performed. Recall in particular the choice of the sequence Xf, see (44). This 
choice was performed on the sequence u^, that had been obtained in the beginning 
of Section 1^ by suitable dilations of the given sequence at P. We focused 

on 2-planes normal to TpE (which was identihed with x (0, 0) by a suitable 

coordinate choice). Fix any ^ large enough and let us read the 2-planes normal to TpE 
back in the picture of uf. we hnd a foliation (whose leaves are 2-disks) of B\^{P). These 
disks are almost normal to £ H B\^{P). The same goes for Q, so for each (. we have a 
foliation (whose leaves are 2-disks) of Bx^{Q). The disks, as ^ gets large, are actually 
staying the same, the only thing happening as ^ grows is that we are only seeing the 
part of the disks inside the smaller balls B\^{P) and B\^{Q). 

For i large enough the disks in B\^{P) and B\^{Q) are small compared to the size 
of B^'^. We then consider a closed 2-dimensional surface, diffeomorphic to a sphere, 
that intersects E exactly at P and Q (and nowehere else - here we are using the fact 
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that S is contained in a Lipschitz graph). This surface should be slightly deformed 
around P and Q so that it is and it contains the 2-disk passing through P and the 
2-disk passing through Q (for a certain i large enough, and thus for all larger i)] this 
is possible by the Lipschitz condition and the differentiability at P and Q. We will 
denote this surface by S and still call it a sphere by abuse of terminology. Without 
loss of generality we can give an orientation to £ H and S and assume that the 
intersection is positive at P and negative at Q. Moreover we can consider a tubular 
neighbourhood of S and we foliate it by similar spheres, each of which contains two 
disks of the disk-foliation. For I large enough all the B\^[P) and B\^{Q) are contained 
in the tubular neighbourhood of the first S. Each 2-disk identifies uniquely a sphere 
(recall that the 2-disks correspond to choices of in Section]^. 

From this (4m — 2)-parameter foliation with sphere-leaves we want to select only 
“good spheres”, meaning the following. Recall the choice of X[ in Section]^ and choose 
the corresponding sequence of spheres from the foliation. Keep in mind that in the 
choice of X^ we had a huge freedom (we only needed to avoid a "^^"^“^-small set of 
values) so the same freedom applies here: we only need to avoid, for each £, a small set 
of elements in the foliation, where smallness is understood in the "^^"^“^-measure with 
respect to (these are bad spheres). Observe moreover that (up to neglecting a 

L^4”*-^-negligeable set of elements in the foliation, other bad spheres) we can ensure 
that for every ^ the restriction ui\s(^ is smooth (by the regularity of the maps see 
Proposition]^ and the fact that there are countably many ui). 

Clearly the spheres Sg tend to the initial S (the one though P and Q). By the 
parametrization of the tubular neighbourhood of S with domain S x (—e,e) we can 
also identify each Sg with S, which will be implicitly understood in (50). 

By the blow-up analysis from Section we have that ug ^ u in away from 
any neighbourhood of S. For every fixed I large enough, we then have two small 2- 
diniensional disks D{P^) and D{Qj) in Sj (the disks from the disk-foliation, we can 
assume that they are small compared to the size of Sj) such that 


Ue\Sj\{D{Pj)UD{Qj)) u\s^\(D{Pj)UD{Qj)) in 

Remark further that, as ^ oo, the disks D{Pj) and D{Qj) are slighly translating 
and shrinking, converging (in any reasonable sense) respectively to P and Q. Moreover 
u is smooth, therefore we have that 


'>^£\se\iD(Pe)uDiQe)) ^ls\{P,Q} in will need uniform convergence only). (50) 

ue\D{Pi) u{Q). 

We have implicitly identified Sg with S in this convergence result. Moreover we 
know from Section that for every £ large enough we can identify inside D{Pg), with 
a finite iterative procedure, a finite number of bubble domains (identified by means of 
a maximal function) in which the energy concentrates: here we can renormalise, i.e. 
reparametrize ug by suitable dilations and stereographic projections. The basic step 
consists of dilating D{Pg) (more and more as £ —)■ oo) and wrapping it on an and 
this reparametrization leads to the first bubble map (pg. For each £ the domain of the 
bubble that we are forming is 5^ \ {small disk} (and this small disks shrinks to a point 
as £ ^ oo). For each £ we replace the small disk D{Pg) with this excised 5^ (domain 
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of (f)i), see the top-left picture in Figure^ On the spherical excised bubble domain 
the (reparametrized) map converges in away from a possibly problematic neck, 
on which we continue the contruction after transforming it conformally into a long 
cylinder (using (46), see hgure top-right). Sometimes we will need to perform this 
basic step on a neck domain (which could have a bubble domain in it)^ Again this 


leads to replacing a small disk in the neck with an excised 5^ attached by another neck 
(Figure [^bottom right). The other possibility that could happen is that on a neck we 
have the formation of a bubble without the need of a dilation, in which case we replace 
a cylinder of hxed size by a 5^ with two disks removec^® compare Figure [l]-bottom 
left picture. After a hnite number of steps of this type we exhaust the bubbles and 
reparametrize on every bubble domain. This leads to the replacement of Si by a new 
domain, that is topologically still a sphere, but looks like a hnite number of spherical 
excised bubble domains S'^^, ..., ^ connected by cyclindrical necks with Si \ D{Pi). 

For D{Qi) we do the same thing and replace it with a connected sum of Sq spheres, 
attached by the neck domains. So Si has been replaced by a new domain Si which is 
topologically still a sphere but, thanks to the renormalization process, the new map 
Ui dehned on it by reparametrization can be passed to the limit in a much stronger 
sense. Observe also that ui is still smooth, since so was ui. Remark that we have 
strong convergence on each spherical bubble domain to a smooth harmonic map, strong 


convergence to u on Si\D(Pi) as in (50) and the contribution to the Dirichlet energy in 
the necks is as small as we like for i large enough (the latter is the content of Theorem 


1 . 2 ). 


To be precise for the forthcoming argument, we need to hll in the missing disks of 
the excised spheres and cylinders, i.e. we want to think of Si as a connected sum of 
spheres (the cylinders also become spheres once we hll in the missing disks): 


Si = Si#^Nij#^Sl, 


with j and k ranging over a hnite set. Here we denote by N the spheres that came 
from the neck domains and by S"^ those that came from the bubble domains. Roughly 
speaking we are thinking of Si as a sum of closed spheres and closed cyclinders, at¬ 
tached by small disks: when we attach them by overlapping the small disks, orientation 
cancellation gives back Si. 

It is not essential to know exactly how we extend the map u on these extra disks. 
The disks can be taken as small as needed, up to chosing (. large enough; by the strong 
convergence on each bubble domain to a smooth harmonic map (a bubble map) and 
by the strong convergence on Si \ U T*^^) to the smooth harmonic map u we can 
dehne u on these small disks by using as little Dirichlet energy as we likf^® say £ in 


total. The area covered in J\f by the image of these small disks is also at most e, since 
the Dirichlet energy controls the area. 

Recall that by Hurewicz theorem is the same as H 2 {N,'L). Consider now the 

map ui : Si ^ AAand remark that (h^)*^*?^]]) is the same integral 2-dimensional cycle 
as (u£)*([[S'£]]). We want to look at the homology class of (u^)*^^^]]). This class can be 


^^Comparin^o Sectionj^this corresponds to the case when we have weak convergence of Ve to Vao. 
^®In Section Im this corresponds to the case when we have strong convergence of Vi to Voo- 
^®We are counting the extra energy brought in by closing the spheres and the cylinders, here each 
disks contributes to the energy twice as there is no cancellation. 
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obtained by adding the homology classes of the spherical snb-domains S'^j^ (the bnbble 
domains) and those of the cyclindrical necks Nij (these are now all closed snrfaces, 
i.e. topologically they are S'^’s bnt we will keep the names to distingnish the bnbble 
domains and the neck domains). 

Claim: the homology class of (nf)*(|I5'£]]) is the snm of the class of n : S —)• AAplns 




Fignre 1: This example illnstrates the passage from to Si at a point of S (say 
P) in the case that there are three babbles arising. In the top-left pictnre we are 
prodncing the first bnbble 4>i. On the part of the sphere above the dotted line we 
have C^’"-convergence to the bnbble map (pi, the neck (annnlns) below the dotted line 
gets conformally transformed ino a long cylinder and we continne on it the babbling 
analysis (top-right pictnre). We are assnming that on this cylinder there is (in the 
centre) convergence to a bnbble (f >2 withont need of reparametrization (this corresponds 
to —)• Voo strongly in with the notations of Section]^. Fo becomes the second 

bnbble (j )2 and we visnalize this by replacing the middle part of the cyclinder with a 
spherical bnbble domain (bottom-left pictnre). Here we see two necks and since there is 
only another bnbble to be tracked down, one of the two necks allows no farther bnbble 
domains (say the top neck). In the bottom neck we assnme that there is a bnbble 
formation bnt it needs a reparametrization (this corresponds to Vi —>■ Vqo weakly in 
with the notations of Section]^. Thns we again replace a small cap with a bnbble 
domain attached by a neck (bottom-right pictnre) and prodnce the third bnbble ^ 3 , 
by repeating the procednre that we nsed for (pi. 
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the classes of Vp (the bubbles originating at P, here p = 1, ...,Sp) minus the classes 
of Wq (the bubbles originating at Q, here q = 1,S'Q-recall that at Q we have that 
D[Q) is oriented negatively): 


Sp Sq 

in ^^ D ^ o 

{ui)4Sij = (u)4>s'l + ^(^^p)4'S' 1 - I- 


p=i 


q=l 


This claim follows by the following facts: 


(i) we have C^’^-convergence on the bubble domains to the bubble maps Vp and Wq 
(thanks to the reparametrization), therefore we can pass to the limit in the ho- 
motopy classes; 


(ii) we have the uniform convergence of Uf\s^ to u\s as in (50), therefore we can pass 
to the limit in the homotopy classes; 


(iii) the energy in the cylindrical neck domains goes to 0 by Theorem 1.2 and a fortiori 
for each neck domain the area of its image goes to 0 as £ —)• oo, therefore the 
(finitely many) necks do not contribute in H 2 {N,'L). Recall that for each non¬ 
zero homology class there exists a mass-minimizer with non-zero mass and there 
exists ^ > 0 such that any integral 2-cycle with mass below ^ must be in the 
trivial class, see |R1| . Here we are also ensuring that the e of extra energy that 
we introduced by closing up the cylinders is small compared to as we said this 
is no problem as long as £ is large enough. 

With the claim in mind, observe further that the map u is smooth on which 

means in particular that u : 5 —?• AAis contractible, so its class is [0] G '7r2(A(). Therefore 


(u44'54 


in 


Sp Sq 


p=l q=l 


(51) 


Moreover the assumption that Ui satisfy ([^ yields as well that (tt£)*([[S'£]]) is in the 
trivial class of H 2 {M, Z), as follows. Denote the bubbles at P with Vp : S'^ ^ AAand the 
bubbles at Q with Wq : S'^ ^ AC, for p ranging in {1, ...S'p} and q ranging in {1 , ...S'q} 
(probably Sp must be the same as Sq but this is not important for this proof). The 
bubbles at P are all holomorphic for the complex structure apl + bpJ + cpK whilst 
the bubbles at Q are all holomorphic for the complex structure uqI + bQj + cqK. 
For any closed 2-form a = aClj -|- bClj + cQk on AC with -|- 6^ -|- 4 = 1 we consider 
(m£)*([[5£]])(q:) = [[5£]](u£Q;) and keep in mind that = dL for some 3-current L 
(since we are in a ball R^™), therefore (u£)*([[S'£]])(a) = L{d{u^a)) = 0 by (|^. This 
is true for all the above choices of a. Unless the bubbles at P and the bubbles at Q 
are holomorphic for the same complex structure (which is the desired conclusion of the 
theorem and in which case we can show even more easily that (u£)*([[(S'£]]) is in the 
trivial class) we can choose a above so that its integration on X^g=i(rr’g)*[[*S'^]] gives 0 
and its integration on X^p=i(4’p)*|I5'^]] is non-zero Since a is closed, its integration 


^°There is an of complex structures on Af and so there is an of complex structures orthogonal 
to agl -h bgj -f CqK and exactly one complex structure that is orthogonal to both apl + bpJ -\- cpK 
and aql -I- 6q J + cqK. 
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gives the same result within a fixed homology class, so we have contradicted (51) 
Therefore the claim is strenghtned to 
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^ ^ o in H2(N,'L) 

p=l q=l 


(52) 


We will conlcude by showing that (52) can only happen if apl + bpJ + cpK = qqI + 
bqJ + cqK, thereby proving the result, Consider the integral 2-current 


Sp 


Cp:=EW-(IS"l). 

p=l 


As a sum of holomorphic spheres it is calibrated by the 2-form apVti -|- bpQj + cp^tx 
(this is a calibration of degree 2 in A/^, i.e. the tangent 2-planes are of the form v A 
{apl+bpJ+cpK)v. Then the standard calibration argument shows that Cp miminizes 
the mass M (i.e. the area counting multiplicities) in its homology class: indeed, let T 
be another 2-current in the same homology class, we then have T — Cp = dL for some 
3-current L in A/" and we can write, using the closedness of apVti + bpVlj + cpQ.K in the 
second equality and Wirtinger’s theorem in the last inequality (equivalently, the fact 
that the comass is 1): 


Ad{Cp) = Cp{apQ,j bpVlj -|- cpVlx) — T{ap^lj bp^lj cp^lx) ^ M{T). 

The amount of information of this simple computation is even greater: equality in the 
last step holds if and only if T as well is calibrated by apQj + bpQj -|- cpQx- We 
conclude terefore that any minimize!' in the homology class of Cp must be calibrated 
by apQ^j -\- bpVLj -|- cpVlx- 

Consider now the integral 2-current 


Cg:= 


Sq 


q=l 


This is calibrated by aq^ti + bqQj + cqVLk and therefore (as before) we have that it is 
a mass-minimizer in its homology class. We concluded in (52) that Cp and Cq are in 
the same homology class. Thus Cq must be calibrated by apQ.i + bpVLj -|- cpi^x, be. 
ap = CLQ, bp = 6q, Cp = Cq. 

Observe also that 0(P) is given by the sum of the energies of the bubbles at P 
by Theorem 1 1.2 [ and for each bubble the Dirichlet energy agrees with the mass of the 
induced current (by conformality). We have seen that M{Cp) = M{Cq) and therefore 
0(P) = ©(Q). 


we assume that on we have smooth maps converging strongly in to ui as k ^ oo, 
then up to neglecting a B^’"~^-negligeable set of further bad spheres in the foliation, this strong 
convergence is valid also for the restrictions to good spheres Se (by the coarea formula). Therefore 
the currents (/fc)*|[S']| converge (weakly as integral currents) to (u^)*|[S']|. By the smoothness of /f the 
maps fl ■. S ^ Af are contractible and so the currents (/lltlSI are in the trivial homology class of 
H2{Af,2) and this passes to the limit by the compactness theorem of |13| . 

^^Remark, as a comparison with the general case of bubbling for stationary harmonic maps, that 
the class of the sum of the bubbles at P is non-trivial, as we noticed in remark [13 there can be no 
cancellation in homology between the bubbles at P (the same happens for Q). 
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As we remarked in the beginning, P and Q can be chosen in an arbitrary manner 
on S n among points where the density of S is 1 and the approximate tangent to 
S exists and agrees with the tangent at C. In particular we have proved that almost all 
points in S n have an approximate tangent for a unique almost complex structure, 
but we still do not know if £ \ S is empty or not. To show that it is empty, and thus 
conlcude the proof of Theorem 1.3 we take two points P and Q' on C with P G S and 
Q' in the open set £ \ S, assuming again that P is a point where TpS exists. Then we 
repeat the argument given in this section. This time the quantization result Theorem 


1.2 is used only at P, while at Q we have no bubbling. Therefore, instead of (52) we 


will conlcude that the sum of the bubbles at P is trivial in H 2 {M, Z), contradicting 
that P was a point of bubbling. We have thus concluded that S n = £ and £ is 
a pseudo holomorphic Lipschitzian graph. 


As remarked in the statement of Theorem |1.3[ a posteriori we have that £ is actually 
smooth. Indeed hrst of all we observe that ^ [£]] is an integral (4m — 2)-cycle that is 
semi-calibrated by the form ( 2 m-i)! {O‘^iPbu}j+CL 0 k)^‘^'^~^\ and thus is has bounded hrst 
variation (we refer to ^ for general semi-calibrations and in particular for those semi¬ 
calibrations obtained by powers of a two-form). Its multiplicity is 1 almost everywhere 
and so we can improve its regularity to by Allard’s theorem [T]. Then we can 
improve it to C°°. 

More directly (even without recalling Allard’s result), we can write £ as the graph 
of / : c ^ m 2 in suitable coordinates. Namely, we can produce a map 

■0 : p4m-2 X p2 gQ for every p G 2-disk V’(P)') is an embedded 

pseudo holomorphic disk with respect to the almost complex structure ai -|- bj -|- ck 
and so that this family of disks foliates P^™'. Such foliation is produced in |3]. We 
can then express £ as the graph of a function / : © 2 ™.-! the 

coordinates induced by 'll;, with /(O) = 0 and V/(0) = 0. Imposing the condition of 
J-invariance on its (a.e. well-dehned) tangent planes, we obtain a perturbation of the 
Cauchy-Riemann equations of the type 


d,f = A]{z, fiz))V,f + P,(£, /(£)), 

where z = {zi,Z 2 m-i) and A^- and Bj are smooth functions that are 0 at ( 0 , 0 ) 
and are small in C^-norm. These two perturbations functions do not have any depen¬ 
dence on V/ in the arguments, therefore this PDE can be treated as a perturbation 
of Cauchy-Riemann by means of standard elliptic theory, and allows bootstrapping to 
C°^. This strategy is an easy case of the one employed in |30| . where a given pseudo 
holomorphic integral 2-cycle in an almost complex manifold of dimension 4 is studied 
and its representation in such coordinates is used to deal with the harder case of high 
multiplicities. 

For the sake of completeness, remark that the smoothness also implies that the 
bubbling conclusions are actually true everywhere on S H P^™', since all points admit 
a tangent plane T^Z- 


7 An application to Fueter sections 

When the domain manifold is 4-dimensional {m = 1) the triholomorphic map equa¬ 
tion (Q corresponds to the classical quaternionic 9-bar equation studied for the hrst 
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time in [15] and thus also called the Fueter equation. In this section we describe a rather 
direct application of our work to the “bubbling analysis for Fueter sections”, treated 
by T. Walpuski |39| . The setting in |39| requires a slightly more general notion than 
the one we have taken in the present paper, namely one wants to deal with a domain 
manifold that carries an almost hyper-Hermitian structure only locally (this can 
be always achieved on a 4-dimensional manifold) but not necessarily globally, since a 
topological obstruction might be present, preventing the global extension of {i,j,k). 
In order to overcome this, one considers a compact bundle x of hyperKahler manifolds 
on Ai with a fixed identification of the unit sphere bundle of self-dual forms on Ai 
with the bundle of hyperKahler spheres of the fibres of x- Moreover one needs to fix a 
connection 1-form A on x- To clarify ideas, observe that in the case of triholomorphic 
maps studied in the present work, we would have the trivial bundle Ai x Af endowed 
with the flat connection and the trivial identification uii — >■ I, Uj — >■ J, — )• K. In the 
bundle one then looks at a (smooth) Fueter section, i.e. u G r(x) that satisfies |39[ 
(1.2) or (B.3)]. We remark, in order to avoid confusion, that the situation studied in 
the main body of [39] would correspond in our situation to a homogeneous triholomor¬ 
phic map (a so-called tangent map), which is why the domain is taken 3-dimensional; 
in Appendix B |39| deals with the 4-dimensional case. As explained in jS^, the study 
of Fueter sections is motivated by Gauge Theory on G 2 and S'pin(7)-manifolds, in re¬ 
lation with codimension four bubbling phenomena. The bundle in question is, in that 
case, a bundle of moduli spaces of anti self-dual instantons on Ai , which can be given 
a hyperKahler structure. 

The presence of the connection changes our equation 0 by adding on the right hand 
side lower order terms (depending on u but not on its derivatives). Indeed, compare 
equation (1.2) or (B.3) in |39j . the term Idui in Q must be replaced with a term of 
the form I{d + A)i, where we are writing {d + A) for the covariant derivative induced 
by the chosen connection. In local coordinates, with summation over repeated indexes, 
this term reads 


/V- a u = Idu 


dx' 


^ dx^ j 


+ 


( f) \ ( 9 ?/“ r) 


d 

dyd' 


Remark that in the last two terms, the first one is just what appeared in (21), so 


the perturbation term in Q is the second term (similar terms appear for (j, J) and 
(fc, K)). We will describe briefly how such terms do not affect the analysis that we have 
performed in order to reach our main results. 

In Section 3.1 we used the triholomorphic equatio n ([I| ) together with the geometric 
structures on the manifolds Ai and AA to infer from (21), in local coordinates on A4, 


the second order PDF (23), which shows the Jacobian structure for A„tt. This was 


the key initial step for the e-regularity result and for the quantization result. In the 
situation of Fueter sections from |39] the perturbation terms do not affect our analysis: 
indeed, using our notations, these terms reflect in an extra L^-term on the right hand 
side of (23) when we compute AgU^, namely the perturbation is 


and the L^-norm of these extra terms is controlled up to a universal constant by the 
Dirichlet energy of u (only first derivatives of u appear and we never take products 
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of derivatives of u). In other words, the PDE structure AgU G /t^ that we produce in 
Proposition on which our later analysis is ultimately based, is preserved with the 
same estimate ||Agu||^i(^) < C" |Vrip (up to changing the constant appearing on the 
r.h.s.). Our analysis can then be carried out in the same fashion in the application under 
consideration. The e-regularity result for Fueter sections was established independently 

In adddition we can give 


m 


and we could also adapt the proof from Section 3.3 


an affirmative answer to the quantization question, which is conjectured to be true but 
left open in 


page 5] and show the analogue of our result Theorem 1.3 


Theorem 7.1. Let ui G r(x) be a sequenee of Fueter seetions, i.e. they satisfy (1.2) 
in m and let us assume that the Dirichlet energies E{ui) are uniformly bounded. Up 
to a subsequence the converge weakly in IT^’^ to a section u (possibly singular) with 
associated rectifiable defect measure (we are using notations as in /391 Thm. 

1.9]). Then there is an a priori -estimate 


< CE{ui), 

for a universal constant C > 0. Moreover for Tl^-a.e. x G S there exist non-constant 
holomorphic spheres (the bubbles at x) Cs : — )■ Xx, where Xx is the fiber on x G M 
and where Sx GN and s = 1,..., Sx, so that the following energy identity holds 

g 

e{x) = Y^E{Cs). 

S=1 

The bubbles at x are a priori all holomorphic for a complex structure that depends on x. 
Whenever S is (locally in an open ball B^) contained in a connected Lipschitz graph L 
(without boundary) and if the limit section u is smooth on then for all x G EClB^ the 
bubbles are holomorphic for a unique complex structure (independent of x) and TiDB 
is a smooth connected curve in B^. 
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